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ABSTRACT 


This final report contains the results of an investi- 
gation to use the conecpts of radiation from accelerating charges 
to develop the radiation characteristics of pulse-excited an- 
tennas in the time domain. The transition to antennas is accom- 
plished through derivation of a radiation equation involving the 
time-derivative of the antenna currents. For sinusoidal time 
variation this equation reproduces well-known results for a 
number of common antennas. For pulse-excited linear antennas, 
the radiation equation produces a result derived previously by 
retarded potential and other methods. Numerical comparison with 
a published frequency-domain analysis of a pulsed dipole shows 
reasonable agreement except for angles close to the end-fire 
direction. Further work on this point is required. For a pulsed 
loop, comparison with a published time-domain analysis using the 
Sommerfeld radiation equation for a small dipole shows that both 
approaches lead to the same final equation. Concerning transient 
aperture antennas an analysis by Cherneusov is described, wherein 
e.\ivalent surface currents replace the aperture fields. The 
Chernousov aperture results are shown to be expressable in terms 
of equivalent surface accelerating charges, and for one-dimension 
to reduce to results derived previously by Cheng and Tseng. It 
is concluded that while the accelerated charge approact provides 
a direct physical explanation of and an analytic basis for study 
ef impulse antenna radiation in the time domain, more work is 
needed to establish its merits relative to other time-domain 


methods. tid 


EVALUATION 


1. Over the past few years there has been an active and increasing 
interest in the behavior of antennas when excited by short 

time duration impulsive like signals and fields. Moreover, 

it is evident that conventional cw concepts and theories are 
inadequate for describing and analyzing impulsive antenna 
performance. The intent of the work reported on herein was 

to examine basic underlying time domain principles and concepts 
useful for understanding and analyzing the basic radiation 
properties of short »ulse antennas. 


2. This report presents the results of a six month research 
effort to develop the concepts of accelerating charges as the 
underlying and basic radiation mechanism for impulse antennas. 

It has been shown that radiation can be conveniently for.ulated 
in the time domain by the time derivative of current or accelerated 
charges over differential elements. Apprupriate expressions 

were derived from basic principles for tinear and small loop 

; antennas which resulted in direct physical explanations of 

z impulsive antenna radiation. Although these same results can 

be obtained by other time domain representations, the accelerated 
charge formulation provides physical interpretation and insight 
into the radiation mechanism for these type antennas. 


3. A promising method for analyzing high gain aperture antennas 
was also identified and discussed. Although not developed fully, 
the method was shown to reduce, giving results to previously 
obtained for a one dimensional case. 


4. Based upon the work performed under this preliminary short 
term effort, it is concluded that the accelerated charge concepts 
provide a valid analytic basis for impulse radiation. However, 
before significant advantages over other approaches can be 
demonstrated, additional effort will have to be expended to 
extend the scope and detail of the work completed herein. Future 
work has been identified and discussed in the report. 
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SECTION 1 


INTRODUCTION 


1.1 GENERAL 

This final technical report covers work performed from 
1 Sept. 1971 to 29 Feb. 1972 under USAF Contract No. F30602-71-C- 
0278/PMRA, entitled Time Domain Impulse Antenna Study, for the 
Rome Air Development Center, Research and Technology Division, 
Air Force Systems Command, Griffiss Air Force Base, New York, by 
Dr. Morris Handelsman, Electrical Engineering Department, College 


of Technology, University of Vermont, Burlington, Vt. 


1.2 PURPOSE OF PROGRAM 

The purpose of this program is to investigate, develop, 
and apply the concept of radiation from accelerating charges to 
antenna systems excited by impulsive signals or fields. Thus the 
purpose is to develop in the time domain, an approach to the 
understanding of the behaviour of impulse-excited antennas, 
through the radiation characteristics of the accelerated charges 


on such antennas. 


1.3 CONTENTS OF THIS REPORT 
In this final report, Section 2 discusses the nature of 
radiation from accelerating charges. Sections 3 and 4% discuss 


the radiation from a pulse~excited linear wire antenna in a 


direction normal to the wire, and in any direction, respectively. 
Section 5 discusses radiation from pulse excited loop antennas. 
Section 6 discusses transient radiation from aperture antennas. 
Section 7 contains conculsions and recommendations. Section 8 
lists references. Seven appendices contain discussions and 
results considered basic and pertinent to the understanding of 
radiation from antennas through the mechanism of accelerated 
charges as follows: I. Lienard-Wiechert Petentials, II, The 
Fields of Moving Cherges, III. Wave Propagation on Tran: mission 
Lines and Linear Antennas, IV. Drift Velocity of Electrons in 
Conductors, V. Derivation of the Radiation Fields of an Electric 
Dipole Using Moving Charges, VI. Some Results of Delta Function 
Integrals, and VII. Illustrations of Use of Equation (74) for 


Some Standard Antennas with Sinusoidal Time-Varying Excitation. 


1.4 PROGRAM ORGANIZATION 

The person who performed this work is Dr. Morris 
Handelsman of the Electrical Engineering Department, College 
of Technology, University of Vermont. Acknowledgements are 
gratefully made to Mr. Hugh C. Maddocks and Mr. Albert E. Ruehli, 
Ph.D. candidates, for discussions cn the Lienard-Wiechert poten- 
tials of the Hertzian dipole, and the frequency spectrum of 


pulses, respectively. 


SECTION 2 


THE NATURE OF RADIATION FROM ACCELERATING CHARGES 


2.1 INTRODUCTION 
' 
—e ’ 
Given the position X (t') of a moving charge | at all 
times ¢ » as shown in Fig. 1, the Fields at an observation point 
POX, t) nay be found through the Lienard-Wiechert potentials. 


These potentials, developed in detail in Appendix I, are in MKS 


units, 
bK,e\= ATE, —\s\. (1) 
Lizt\= aterm on 
RGIS aa LS Saat 
where 
S=r- —e ~y—-¥. 8 =v (\-Acose) (2) 
velocity vc. = ar / at’ (4) 


_— 


¥s) alc (5) 


and E | + means that the quantity inside the brackets is to be 
ve 


iH 


evaluated at the retarded time 


tot- Vie (8) 
t =/ 
Thus t is the time at which a signal is emitted at X so as to 


arrive at x at time ¢ 


(t'), 


FIG. 1 Charge q at retarded position X 
4 


2.2 THE FIELDS OF MOVING CHARGES 


: The fields B and E are given by 

, aA 
Ex si 2 (8) 

at 


This results in the following equation for the fields (for de- 


tails see Appendix II): 


= A (¥-8 1) C-.3? 
e at (¥-8 r )C\- 9? ) 
+ oreo tatesoesli) 
8=ahe LF (\ 8”) aeae Nix S0xa]} 
ref, 
(10) 


Here & = acceleration of the charge -aw/dt 


For a stationary charge, Eq. (3) shows that OS=V , 


Cal 


while B and @ are zero. The E field then becomes the static 


Cuulomb field 
A 


ee a 
CS ate 


A ~~ 
where Gr = unit vector in ¥ cirection, while @=0. Fora 


Sneed 


charge moving with uniform velocity without acceleration, p= 
constant and @=0. The E and B fields vary as ifr? » and are 


non-radiation fields (quasistatic or induction fieids). 


2.3 THE RADIATION FIELDS 

When the charge has acceleration a, the last terms in 
Eqs. (9,10) are the radiation fields, vary as 1/¥, and depend 
upon a. These radiation fields, denoted by subscript a for 


acceleration, are 


eo [ fastierna) (12) 


a 2 3 
4NE.C S cet 
B = reife LG ¥xire(G-Snxalh 7 (13) 
© 4ié, meee 3\ 
tet 
Inspection of Eqs. (12,13) shows that 
9 f —. 
8, = 4 Ay xX e. (14) 


Hence 
W=atS8= yet, 4 ce (15) 
Wa = un = Eq [tl Axe, 
a oe 8 ; 
That is, Ea, Ha, and &,(the retarded position unit vector) are 


mutually perpendicular, and the ratio of \EQ| to (Ha 


Ea j Hale J Mo 1€. ohms cte7 


which are all familiar relationships to the antenna engineer. 


2.4% LINEAR ANTENNAS 
For a linear antenna extending along the Z-axis, as 
shown in Fig. 2, velocity v and acceleration a are colinear. 
Then, in Eq. (12), 
(¥-Sr) ka = TxA (a7) 
Hence, omitting the U JV potation, 


a mh Lexa) 
AT Eo c*s3 


gE. Ses, (eee Geka era) _ (18) 
4TE,c* FCI- Bcose) 3 


fn 


i 


Han dex G (19) 
J Wel, 
—_ A 
Using axa, a 
A = 
0. 4(A Ka) = AMS ae (20) 


A 
where ae, is the unit vector in the increasing 9 direction, 


using spherical coordinates. Then 


a, SiIN® (21) 
4 & c* ¥ U- 8 @se)\> 


He - Ee / [acTee (22) 


= al 


Linear antenna, 


2 


FIG. 


eon 
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aera ee 
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The radiation patterns which result from Eqs. (21,22) 
are shown in Fig. 3. For low velocity particles, which is the 
usual case for charges (electrons) moving in conductors (see 
Appendix IV) for detailed discussion on the drift velocity of 


electrons in conductors), Pae\ , and 


- 


me = re a ( Beet) (23) 
41E, C70 
ret 


This is the "figure eight" pattern, shown in Fig. 3(a), which 
varies as sin 8 » again familiar to the antenna engineer. For 
relativistic speeds f-»>| , Eq. (21) shows that there is a 
“forwari" bunching, as indicated in Fig. 3(b). The direction of 
maximum intensity is now tipped forward. 


The H field is 


- 


Hy — |—Yasme (Bee (24) 


cv 
AT tet 


To illustrate an application of Eq. (24) consider the 
radiation field of an electric dipcle or current element of len- 
gth \ » extending along the Z-a.is, with a current distribution 
along \y which is constant except at the ends where the current is 
zero. oce Fig. 4. Let the current vary sinusoidally in time. 
This case is examined in more detail in Appendix V. The sinu- 
soidal case is examined initially principally because the H 


fields for this case, using conventional theory, are well-known. 


aq 


(a) Velocity vege 


(b) Velocity v approaching c 


FIG. 3 Radiation patterns of accelerating charges. 


to PCT) 


FIG. Y Electric dipole or current element. 


The H field for h small compared to wavelength r » is (Ref. 20, 
p.- 498; Ref. 21, p. 93) st 
: . — JRC Jue A 
i = _\s Tolan Md eT Co 4% (25) 
TY 
where K=wW(c=2T/) , and the current by definition, is 
( 
ae \wt a 
T=r,e) as (26) 
The charge dq in a differential length AZ at ( o, o,2 ) on the 
radiat ris 
dq = PA d?2 (27) 
where P is the moving (electron) charge density in coulombs per 
cubic meter, and A is the appropriate cross-sectional area. 
( 
This charge moves through distance dz in time dt » so the 


current is 

_ A f . 

T - 24. a, = PA oz A, = PAV (28) 

ot t! 
— ‘ A 
where VU =(02/dt') is the velocity of the charge. The current 
density (amperes per square meter) J is 
J aL Ser oe 

Assuning P is essentially constant in o7c2ch (see Appendix ill! 


— 


for discussion), then the acceleration @ is 


ee | oF (30) 


From Eq. (26), 


es _ of 
) joT ci (31) 


dt 
Replacing q in Eq. (24) by the above dq » then the term (qa ) 


in Eq. (24) becomes, using Eqs. (27,30,31), 


qa= eae = = jot, azo” (32) 


( 
The retarded value of (4&) is found by using Eq. (6) fort , 
so that 


Since dq is at distance 8 from P at time +' » the T in the de- 


nominator of Eq. (24) becomes R. Thus dWe due to dq is 
AUG t= joTactt gies eS sne! de 
41 


and 


LY MN — 
He t= ert | CMR sino’ yo 35) 
R 


Te 
to) 


For the radiation field (far-field distances) 


t 
620 
(36) 


RY-2@50 


2419) Ye 


Using Wee) » the integral in Eq. (35) reduces to 


Wr 42 P 

j _@ SER sine! ya - ea 3\n0 (37) 
R 

fe) 


Hence the field calculated from Eq. (24) becomes 
s ‘ —Jke jot 
Wy N= VSEoh_gmo6 e (38) 
we 


which agrees with Eq. (25). 


2.5 HOW RADIATION WAVES ARE PRODUCED BY AN ACCELERATING CHARGE 

It is possible to derive the radiation equation Eq. 
(23) using a graphical construction and simple analysis (Ref. 
23, pp. 60-63 or Ref. 24, pp. 334-340). This derivation may be 
of help in understanding how the radiation field of an accelerat- 
ing charge arises, and is presented here. The charge 9 is taken 
to move on a straight line. Three preliminary points are made 
first, in the paragraphs below. 

First it is noted that the term (F- BY ) appears twice 
in Eq. (9) and therefore deserves interpretation. Remembering 
that ¥ is the retarded position of the charge, let Ge be the 
"virtual" present position (at time t) of q, i.e., the position 


it would occupy at time € if the acceleration Gis zero. In 


Fig. S the time of propagation from the retarded position Q, 


*3 


Mek ite Pity, MA” 


FiG. 5 Virtual position of charge moving at 
constant velocity along straight line. 


GI 
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where the charge is at time t-Vic to the observation point P is 
simply 
t=atic (39) 
During this time t » the charge is moving with constant velocity 
VizBc, and zero acceleration, along the line, so it reaches the 
point Q at the same time that the energy which left from ie 
reaches P, 
Hence 
Q.Q,= Tt= Hr (40) 
It is seen that 
V2 yer (41) 
Thus the term (Y-SY) is the virtual present position. 
Second, the quantity S in Eq. (3) can be identified on 
Fig. 5. Since 
~- 8 =( Bl) we (42) 
then (#8 ) is the projection of the vector Mr on ¥ , ice., 
the segment QT of the line ¥ = Q,P. Then from (Eq. (3), it 
follows that 
S=r-"- B= TR ae 
Further, S can be written in terms of Vy and Oy » the virtual 
position and angle, respectively, of the charge at q, . Thus in 


the right triangle Q ve, 


2 2 
S*= Yay — A rsin*@ 


But the perpendicular H from P to the line of motion is 


W= Tome = Ty aN by (44) 
Hence 
2 2 
(a9 = ty SI Oy (45) 
and 
f2 
S ‘as Vy (\-." sin“, ) (46) 


Third, the first term in Eq. (9), designated as E 


can now be written, using (Eqs. (41,46), as 


= ig -g2 
Ee =| Ate os - - (47) 
| 4Teo Gy | | (l- 975im79,)?? 


where One is the unit vector along Vy » pointing from the pre- 
sent position of Q to PF. It is seen that Ey» which is the 
total field of a charge moving at constant velocity, since accel- 
eration & is zero, is identical to the Coulomb field of a sta- 
tionary charge, given by the quantity in the first brackets of 
Eq. (47), modified by the quantity in the second brackets, which 
is a function of the velocity. Eq. (47) agrees with Ref. 25 

(p. 254). For low velocity particles, where fiée| , the quantity 
in the second brackets is essentially unity. Hence the impor- 


tant result: The field of a charge moving at a low constant 


velocity is essentially the ordinary Coulomb field associated 


ee ED 


with a stationary charge. Thus, as put in Ref. 24 (p. 337), for 
a charge moving with constant velocity, the lines of force di- 
verging radially from the charge move with the charge, acting 
as if they were rigid wires attached to the charge. 

Now the effects of acceleration can be ascertained. 
Following Ref. 24 (pp. 334-340), let point charge 4 move along 
a line as shown in Fig. 6, with constant velocity W2é4c, arriv- 
ing at point 0 at time aan . During a time interval At, let q 
be accelerated to velocity V+AU, reaching point O, at time 
Lorot » after which it moves at the new constant velocity. As 
established above, the E lines move with q as long as it does 
not accelerate. The effects of any disturbance such as the 
short period of acceleration are assumed to propagate outward at 
the speed of light c , which is reasonable. Draw a sphere S, 
around 0 of radius 

You c(t-te) (48) 

where tot, . Any signal traveling at velocity C reaching any 
point in the region outside S, at time t must have left Q ata 
time before x, . Thus the E lines outside S$, are the same as 
those of the charge moving with the original constant velocity 
a. Therefore, by the third point established above, the E 


lines outside Sg are essentially the same as these ‘+ charge 


FIG. 


6 


Radiation due to an accelerated charge. 


a PG 


at the point Qy, which is where the point where 4 would be at 
time t if it had not been accelerated. Hence the E line BC 
outside S, is shown diverging radially from Q, . 

Next draw a sphere S, centered at 0, with radius YQ - 
cat, Any signal traveling at velocity C reaching any point 
inside S, must have left q after it reached 0, at time totat, 

at which time it is moving with velocity W+OU. Hence the 
field inside S, is the same as that of a charge moving at con- 
stant velocity Trav, and therefore diverges aS the line Qa 
from Q, ; which point is actually reached by q at time Tt . 

The thin shell of thickness Cov is the region where 
the acceleration effects must put a "kink" in the E line Q, ABC 
as shown in Fig. 6. This kink travels along the E line at veloc-~ 


ity C. Then, 


- 4 (49) 
ae at 
60, 3 Wot (50) 


Since W/lczc\ , then 00, 22 Ct. Hence the shell is approxi- 
mately due to two concentric spheres, with coincident centers at 
0. To an approximation which negi.cts terms of the order of 
(asl) as well as (ar? (ce? ) » the situation is as shown in Fig. 6, 
with Q,4 % No (Ref. 24, p. 340) in the limit of very small 


(lc), 


The E field segment AB, which represents the accelera- 
tion effects, is resolved into two components: 


Eo _ AD Qo Gi Sine 


=. = (51) 
Ey DB cot 
But 
Q,p Q, = 08,- OQ, % (4-teh(wrar)-(t-t yw 
— (t-t.\ov 
Qa. % Nyod5 _ Yoaot (52) 
: c Cc 
Hence 
Ea. ~ _Yoasine (53) 
E, c+ 


Continuity of the E, component withthe usual radial field de- 


mands that 


E,= eet re (54) 


Hence 


E,= ne - (55) 
48, C2, 


which agrees with Eq. (23) for the radiation field of an accel- 


erating charge with velocity much less than c. 
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SECTION 3 
RESPONSE OF A LINEAR WIRE ANTENNA TO PULSE EXCITATION 


AND RADIATION IN THE BROADSIDE DIRECTION 


3.1 TRAVELING PULSE ON WIRE ANTENNA; THE SOURCE OF RADIATION 
A current pulse traveling on a linear wire antenna is 
shown in Fig. 7. The current 1(Z,t) satisfies a wave equation, 


assuming no attenuation for simplicity, given by (see Appendix 


IIT) “ 


ae eae se OTe (56) 
oZg* <c*  or* 
The solution to Eq. (56), for a wave traveling in the +Z direc- 
tion is the well-known equation 
T(42= T(+- ze) (37) 
For ordinary conductors, it may be assumed that the 
wave velocity is essentially the velocity of lipht c (see Appen- 
dix III). Skin effects will produce some round-off of sharp 
wave-fronts. In the central part of the pulse shown in Fig. 7, 


where the current is constant, no radiation is produced (a d.c. 


current does not radiate). This foliows also from 


~— 


TL = TAHPVA = constent = TI, (58) 
Thus velocity v is constant, acceleration a is zero, and there 
is no radiation from that portion of the pulse where the current 


is constant. 


T(z) 


FIG. 7 Traveling current pulse on straight wire. 


J=constant= i. 


fri FIG. 8 Leading edge of traveling trapezoidal pulse 
: shown at time 7 . 


The radiation arises solely from the leading and trail- 
ing edges of the pulse, where the current (density) is changing 
with time, which results in acceleration of charges. To sim- 
plify the discussion of this phenomenon, consider an idealized 
trapezoidal pulse with a very short rise time = tp seconds at 
its leading edge, as shown in Fig. 8. This figure shows the 
pulse at time t, (i.e., a "shapshot of J "frozen" at time 4: ). 


In time 


terdz/c (59) 


J at point O changes from zero to J; . Thus the partial deriva- 
tive of J with respect to time is 


oT _ Werte, % \- I(t 2) . T-0 . Te (60) 


ot” te te te 


In Eq. (60), Z2 is fixed at its value 2, at point P. The partial 


derivative of J with respect to distance Z is 


ay ves AV G9 =) Ct -O2) _ O-3o i -Se_ (61) 


oz 2 S22  ~ a2 


From Eqs. (59,61), 


Differentiating the equation Ju PV with respect to time results 


in 


oT pw iis OP (63) 


The continuity equation, for this one-dimensional case reduces 


to 


- 2£— 20-3 = oq (64) 


Combining Eqs. (60-64 incl.), the result is 


ay or ws 
- + (65) 
te c Of Cte 


Since v, the drift velocity of charges (electrons) in a conduct- 
or is many orders of magnitude less than c (see Appendix IV), 


Eq. (65) reduces to 


ar. & j- S)y He 
or ogee eal c =e (66) 


Thus, alchough aPlat can approach large values for very smal] 

te , the term WOMdt turns out to be smaller than Povlat by a 

factor equal to the ratio of the drift velocity to c, which is 

considered to be on the order of io” or less (see Appendix IV). 
The above results for the trapezoidal pulse are a 


special case of the more general result discussed in Appendix 


Pan /p ov oy (67) 


MO oe 
v 


or‘ ot’ 

Hence, using Eq. (27), for charge q in length dZ, 

[aa] sO TA aane vee (68) 

Tet ot 

which is the same as Eq. (32). Thus the radiation arises from 
the time derivative of the current, i.e., from the edges of a 
flat-topped traveling pulse, as previously asserted. Insertion 
of Eq. (68) into Eqs. (23,24), and integration over the antenna 
length variable Z, taking into account the retardation factor, 
allows calculation of the radiation fields. Further discussion 


of this is given in Section 4 (e.g., see Eq. (74)). 


3.2 RADIATION NORMAL TO A WIRE ANTENNA 

To illustrate the use of some of the concepts concern- 
ing radiation which have been presented, they will be applied to 
a well-known paper concerning the calculated and experimental 
response of a thin cylindrical antenna to pulse excitation 
(Ref. 26). In Ref. 26 the radiated field of a monopole at B=90 
is calculated by Fourier transform methods, and confirmed by 
measurements. The monopole, of height \ over a conducting 


ground plane, is driven by a 50 ohm coaxial line, as snown in 


Fig. 9. The antenna length-to-radius ratio (h/Q@ ) is given as 
904, Several pulse lengths are used with space lengths CY vary- 
ing from 0.2 h/c to 2.0 h/c. The exact shapes of the pulses 
are specified in Ref. 26. In this report these shapes will be 
taken as almest rectangular, with short rise end decay times, 
which greatly simplifies the discussion. Further, it will be 
assumed that tne input surge impedance of the monopole at its 
base is approximately 300 ohms. This prod. s a current reflec-~- 
tion coefficient of 0.72 on a 50 ohm line, which, as shown 
below, gives good agreement with the published results. 

The calculated radiated E field (Ref. 26, Fig. 4, 
top graph) for T% 0.2 h/c is sketched approximately in Fig. 10. 
The amplitudes have been normalized so that the first pulse has 
unity amplitude. The shape of this graph will now be explained. 

As the leading edge of the pulse emerges from the 
coaxial line onto the monepoli, there is a reflection back into 
the coaxial line, which is of 10 concern here, since that line 
is terminated in 50 chms. As explained in Section 3.1, the 
leading edge of the pulse traveling on the monopole produces 
radiation due to accelerated electrons from the instant it 
energes at local time tiso » Which reaches PCR, 9° ) at time 


t= t+ Ric =Ric. As the leading edge proceeds to a distance 


Yr 
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FIG. 9 Monopole over ground plane. 
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FIG. 11 Radiation from leading edge of pulse. 


Z=ct', which is less than the space length of the pulse, radi- 
ation continues to be emitted from the leading edge, arriving at 
P at time t = t'+Rle, See Fig. 11. Hence the radiation emitted 
in a period of time +! is received over the same period of time 
at P. Neglecting line losses and the attenuation of the pulse 
due to radiation, this radiation continues to arrive at P from 
the leading edge, until the trailing edge of the pulse emerges 
at time t'= T from the coaxial line. At this instant, radi- 
ation being emitted by deaccelerating electrons at the trailing 
edge, with radiation opposite to that of the leading edge, 
leaves to arrive at P at time T+RI¢C, where it cancels the 
radiation from the leading edge. Hence the radiation at P con- 
sists of a pulse of normalized amplitude equal to unity for a 
time interval equal to YT, and zero amplitude afterwards, during 
the time that the entire pulse is traveling up the monopole, and 
before the leading edge reaches the tip of the monopole. This 
state of affairs is illustrated in Fig. 12, where the common 
time delay of RIC has been dropped, and the traveling pulse has 
a normalized amplitude of +1. At the leading edge, charge is 
accelerated from zero velocity to unity (normalized) velocity, 
while at the trailing edge, charge is deaccelerated from unity 


to zero velocity. This is illustrated in Fig. 13, which shows 
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(a) Pulse on monopole (b) Received waveform 


FIG. 12 Initial radiated field when entire 
pulse is on monopole. 
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FIG. 13 Cancellation of radiation from leading 
and trailing edges. 


the resultant cancellation between the two edges. Due to the 
ground piane, all of the above currents and charges have images. 
However, for @=%0 , the images merely duplicate and double the 
radiation produced by the real monopole, so image effects need 
not be considered separately. 

When the leading edge reaches the monopole tip at2ésh ; 
and time t =hle » the current goes to zero at the tip (approxi- 
mately, since there can be c.irrent flow over the end cap, unless 
the end is sharpened to a needle-point (see Ref. 27)). Thus the 
pulse breaks up into an incident and a reflected pulse, which 
must appear as shown in Fig. 14 at some instant after the lead- 
ing edge has reached the tip, but before the trailing edge has 
reached the tip. The combined space lengths of both pulses is 
CY. The radiation from the leading edge of the reflected pulse 
and the trailing edge of the incident pulse reinforce each other; 
hence the amplitude of the radiation field during the reflection 
period of duration T is approximately -2 (normalized). That is, 
the field amplitudeX%-2 in the time period ets R47 » omitting 
Ric. This reinforcement is illustrated in Fig. 15. The trail- 
ing edge of the incident pulse deaccelerates charge from unity to 
zero velocity, while the leading edge of the reflected pulse ac- 


celerates charge from zero to negative unity velocity. In fact, 


FIG. 14 Pulse after reflection at tip. 
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FIG. 15 Reinforcement of radiation from edges 
during reflection from tip. 


the reflection from the end is probably not perfect, and a cur- 
rent reflection coefficient % -0.9 leads to a field amplitude 
=a .\.9 which gives better agreement with the calculated field 
amplitude 2=-\.93. 

Summarizing the above, there is a burst of radiation 
as the pulse emerges from the coaxial line, which ceases and 
remains zero while the entire pulse length is traveling on the 
line. This is followed by another burst of radiation when the 
leading edge reaches the monopole tip, which ceases when the 
trailing edge completes its reflection from the tip. After the 
pulse has completed its first reflection from the tip, there 
results a negative current pulse of amplitude -O.D traveling down 
the monopole, as shown in Fig. 16. Again, as previously describ- 
ed, the radiation from leading and trailing edges cancel, leading 
to zero radiation. Hence the radiated field appears as shown in 
Fig. 17. 

When the reflected pulse shown in Fig. 16 reaches the 
base of the antenna, it is leaving a line with surge impedance of 
(approximately) 300 ohms and encountering a coaxial line of char- 
acteristic impedance =5SOohms. Assuming for simplicity that the 
surge impedance is largely resistive (there can be an apprecia- 


ble reactive component; see Ref. (28), Figs. 5 and 6), the 
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FIG. 17 Radiation field after first reflection. 
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current reflection coefficient is 

Be & =, BOSS Lo y2 (68) 

Ret RL Boo t5o 

The situation during this base reflection period is shown in 
Fig. 18. The magnitude of the reflected current pulse is (-0.9) 
(0.72) = -0.65. The trailing edge of the incident pulse deaccel- 
erates charge from -0.9 to zero amplitude, while the leading edge 
of the reflected pulse accelerates charge from zero to -0.65 
amplitude. The radiation field subsequent to the reflection of 
the entire pulse from the base therefore appears as illustrated 
in Fig. 19. 

The reflected pulse of amplitude -0.6§ travels back up 
the antenna, and the previously described reflection phenomenon 
at the monopole tip reoccurs. Assume Coe O.D at the tip; the 
reflected wave amplitude is 0.59, and the situation during this 
reflection is shown in Fig. 20. The lagging edge of the incident 
pulse deaccelerates charge from-a68 to zero velocity, while 


the leading edge of the reflected pulse accelerates charpe from 


zero to 0.59 velocity. The resultant radiation is then -(-0.65) 
+ (9.59) = 1.24. The radiation waveform now appears a: shown in 
Fig. 21. 


Carrying the process through one more step, the reflec- 


ted pulse of amplitude +0.59 returns to the base, where using 
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FIG. 18 Incident and reflected waves at bas. 


FIG. 19 Radiation field subsequent to reflection 
from base. 
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FIG. 20 Second reflection from tip. 
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e = 0.72, another reflected wave of amplitude 0.43 starts 
back up the antenna. The situation during this second reflection 
at the bas# «ss shown in Fig. 22. The radiation waveform is shown 
in Fig. 23. 

Comparison between the radiation waveforms of Fig. 23 
and the calculated result (Ref, 26) shown in Fig. 10, shows a 
fairly good agreement. %f course, the surge impedance and re- 
flection coefficients have been selected, using previously pub- 
lished data, so as to obtain good agreement. Nevertheless, it 
is evident that the principal radiation characteristics of the 
pulsed monopole antenna, at least for B=% , can be explained 
with the aid of the radiation mechanisms described above and some 


simple transmission line calculations. 
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FIG. 22 Second reflection from base. 
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FIG. 23 Radiated waveform subsequent to 
second reflection from base. 
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SECTION 4 


RADIATION OF A PULSE-EXCITED LINEAR WIRE ANTENNA IN ANY DIRECTION 


%.1 DERIVATION OF EQUATIONS 

In this section, the equations are derived for the ra- 
diation fields of a pulse-excited, linear, thin-wire antenna, in 
any general direction at angle @ with the antenna, using the ac- 
celerated charge concepts previously developed. See Fig. 24% 
which shows a standard center-fed linear antenna. However, the 
following discussion also pertains to any type of linear wire 
antenna, including a monopole over a ground plane. These equa- 
tions agree with those given by Manneback (Ref. 29, Eq. 12), 
Schelkunoff (Ref. 30, pp. 102-109) and Ross et al (Refs. 31,32). 
For an infinitely thin wire antenna, Refs. 29-32 inclusive give 
the radiation field of a general current wave I (Z,t) traveling 
undispersed at light velocity c in the positive Z direction, cri- 


ginating at origin 0, as 


Hp = (t-te) \teose (70) 


47r Sine 


Eg = lle Hp (71) 


An alternative form of Eq. (70) also appears, using the trigono- 


metric identity 


(cose a Sine (72) 
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FIG. 24 Linear antenna geometry, 
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which is 


H = TC +-vic\ SMe (73) 
nn Ts \-cos 


From the accelerated charge viewpoint, the radiation 
field due to an accelerated differential charge q of length dZ 
is found by combining Eqs. (24) and (68) to give 

Att — __ Sine OL CZ, t-R/c) d2 (74) 

4m7¢cR 3 ot 

In Eq. (74), R appears in place of r as used in Eqs. (24) and 
(68). In Eqs. (24,68) the distance from radiating element dZ 
to field point P was designated as r, whereas in order to conform 
to Refs. 29-32 inclusive, r will be used to designate the radial 
distance from origin 0 to P, and R will designate the distance 
from dZ to P. 

The radiation field given by Eq. (74%) is also that of 
a Hertzian dipole, as given by Sommerfeld (Ref. 19), in the form 
shown as Eq. (V-2) in Appendix V, where the dipole length 1 = dZ, 


and 


az 


(75) 


oo 2 _ a 
6(2,t-RIcV= 25 (4(2, + Ric)Az) = BTR c) 


An equation equivalent to Eq. (74) also appears in Manneback 
(Ref. 29, Eq. (4)), where it is attributed to Hertz. In Appendix 


VII, it is demonstrated that the use of Eq. (74) for a number of 


standard antennas with a sinusoidally-time varying currents, such 
as a dipole and a traveling-wave antenna, leads to agreement with 
well-established results. In addition, in Appendix VII it is 
also shown that for arbitrary time variation, dH obtained from a 
vector potential A formulation is the same as given by Eq. (74). 

Consider the field of a current wave traveling at ve- 
locity c in the positive Z direction on the upper half O£ z<\ 
of the antenna shown in Fig. 24. The field due to waves on the 
lower half -\ £2£0 of the antenna is calculated similarly as 
shown later. From Eq. (III-13), the traveling current wave is 
given by 

D(z t)= TU-2/e) e787 

In order to postpone, for the time being, inclusion of the ef- 
fects of terminal conditions at the end of the antenna 2 = h, 
(i.e., a reflected wave moving in the negative Z direction), let 
I at any particular time extend from 2 = 0 to Z= L, where LGC h, 
as shown in Fig. 24. The field at P is then found by integrating 


Eq (74) from 2 = 0 to Z = L, 
L 


Wo eu 308 0 T(2,t-RIc) dz (77) 
e- ane J Rot 


For point P sufficiently distant, the variable R in 


Eq. (77) may be replaced by the usual tar-field approximation 


RX ¥- 2@9 © Te) 


Eq. (78) is based upon the assumption that at any time the largest 


spatial extent Z = L of the current on the antenna is small com- 
pared to R and r, and Eq. (77) is based upon the assumption that 
the current wave I has not yet reached the antenna end Z = h. 
Then, and only then, R may be replaced with negligible error by 
r in the denominator of the integrand i. Eq. (77), and sin ® may 
be regarded as constant, since R and r approach parallelism when 


Eq. (78) holds trne. Then Eq. (77) becomes 


L 
~ sane. DB _T(2t-Ric) d2 (79) 
4 4ire ) ot 


The function I (%,t-R/c) in Eq. (79), can be rewritten, 


using Eqs. (76,78), as 


Zz 
T (at RST 2-B\=Te- F- Fld) 


(80) 

Thus Eq. (79) becomes 

: aoe 
WH, — .2i6.} 0 T(t--lh-+ \- @96\\ d? 
¢ 47 ¢ | ot Cc Cc ( J Goan 
Now 

/ 
OT ({t- L.--2 (ese))= TS t- $-2(\-ase))=— 
ot (82) 


/ 
where I denotes the derivative of I (t- LL ~ Z(\-cose) ) with 
respect to the argument ( t- L- Elune). 
Similarly 
( 
2TK- & - Z(-ws))— ~ me T (83) 
Cc Cc c 
oz 
From Eqs. (82,83), 
a t = 

aT (e-L- 2-9) = = S_ Atlt-F- Fwey) 

at \- 856 Oz (84) 
Eq. (84) is typical of the relationship between the t and Z 
partial derivatives for traveling waves. Substitution of Eq. 


(84) into Eq. (81) gives 


L 
—~_ se _ | oO T(t- L_ 2 
4 = pe alee C ae (\-@56)) 42 (85) 


The partial derivatives of I with respect to either 2 
or t cannot have infinite values, in the actual case. The slopes 
of any real wave I traveling on a wire may be very large, but 
must be finite, even for a hypothetical square pulse applied to 
the input terminals. This is discussed in more detail in Appen- 
dix III, subsection 3, entitled "Wave Travel of Step Functions 
and Pulses on Lossy Lines." However, the employment of step and 
delta waveforms in concept, or as models of actual waveforms is 


so useful and important that their radiation properties must be 


Le - 46 - 


analyzed. This is done in Section 4.2. At this point, assuming 
that I has no infinite slopes, the integral in Eq. (85) reduces 


simply to 


Welt i) _lra- Wley-U(k-vle-b (eesevie| 


oe (86a) 
F(6) = SiO s,s _ tas © (86b) 


\-we ~ se 


The retardation time r/e identifies the first term in 
the brackets of Eq. (86) as the current that was at the origin 


0 at time r/c ago. The retardation time 


Es fe - J-L@e (87) 
ere (rene) - Cc 


of the second term is the sum of the time L/c for the leading 
edge of the current wave I to travel along the antenna from the 
origin to 2 = L, plus the time (v-Leese }(C for propagation of 
the field from Z = L to point P. This identifies the second term 
as the current that was at the origin 0 at time ( UL Giness\ | 
ago, which is the current at the leading edge of I. But for 


the use of non-infinite derivatives of I, the leading-edge cur- 


rent must be ze:%, as illustrated in F 25. Hence Eq, (86) 
nae -~ uy - 
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FIG. 25 Leading edge detail. 
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FIG. 26 Finite-slope function Tcat') 
at fixed time t! 
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reduces to 

H(t") = Lt-tlc} _\twse (98) 

Oo 40r 2nd 

Eq. (88) is the same as Eq. (70) or Eq. (73). The 
field appears to be radiated in spherical waves from the origin. 
Thus, although the radiation emanates from all parts of the 
traveling wave where there is acceleration of charge (i.e., a 
non-zero time derivative of current), the integration of this 
derivative along the wave produces the interesting result that 
the far-field radiation acts as if it emanated only from the 
point of initial excitation, the origin. Similar results, but 
for differing reasons appear in Refs. (29-32 inclusive). 

Eq. (88) can be derived in an alternative manner, 
which is shorter and instructive. In Eq. (85) make the follow- 


ing substitutions: 


{'= +- r/c (89a) 
Zz = Z (\-@se) (83b) 
L' = L (ese) (89c) 
Then "y 
/ 
ay yy— 2sime 1 ( a Tt )da (90) 
Hy (tedor) = aan sar | oe C 


/ : : : 
Let I(t +2 ) be a finite-slope function, as shown in Fig. 26. 


Replace tt integral in Eq. (90), using equi-valued intervals 


ve - ug - 


pays. 
aS Bete aes See 
RT ptntetreeren cine a se sae er 


ae given by 
A?z=uiL/mM (91) 


by summation shown below: 
( 


L (ley t coe ‘ 
a 2 —-FS 2 It-zilc)\o 
{Bret Zic\d2' = F 3, LU 2/c\o2 


A mM 


- L- To) , Y2-T, ace oek: O- Tet a2! 
a2" a2! a2! 


Hence 


L 
| 2TH-2ie\dzZ --Tot' (42) 
d2' 


; ‘ : 5 os : 
Since I (0,t ) is the current at the origin, i.e., 


' m\u 93 
T (22¢,t') = Tlt- tle) (93) 
then Eq. (90) reduces to Eq. (88), i.e , 
Ha (t ‘\- Tht-tic) lteve (sas 
: 470 Sine 
4,2 RADIATION DUE TO STEP-FUNCTION EXCITATION 
The above method, based upon Eq. (90), also allows a 
Simple derivation of the radiation field due to current waves 
with infinite derivatives, such as step-functions, rectanpular 
pulses, or impulses. The transformation from the time deriva- 


tive of Eq. (81) to the space derivative of Eq. (85) continues 


to hold, even for these cases, as can be shown by a limiting 
process going from non-finite slopes to infinite slopes. Thus 
Eq. (90), based upon Eq. (85), continues to hold. 
Consider the traveling-wave step function 
TC2,t\= TL Uit- 2lc) yee 
where U(x) is the standard step function given by U=l for x?0 
and U=0 for x40. Then Eq. (90) yannes 


Us = -Sine | [ afrovue- Z'Ic))d2' (96) 


\-co6 9) «64: 


The derivative of the step function inside the inte- 
gral in Eq. (96) yields a delta function, which can be inve- 
grated. However, the integral can be evaluated alternatively by 
first allowing the traveling-wave to have a steep, but finite 
slope at its leading cdge, performing the integration, and pass- 
ing to the limit of an infinite slope. See Fig. 27. This mathe- 
matical process, has real physical significance, since it shows 
how a hypothetical step function can be approximated as closely 


as desired by an increasinply steep but finite-slope wavefront. 


Hence 


7 
"he 


2 (LOU Zc hz eats az 2 Tt 42 57) 
} 32 


T(z. t') 


I 
\ 
true step | 4 Pup Eros 
To) infinite Sorel \ 
lim a2-oj i \ 
' \ 
\ 
6) L'-a2" dee Lo 2 
FIG. 27 Step function as a limiting process. 
—_—? 
—h ws To pout ee 
qc 


ee an at emer aan 


FIG. 28 


Geometry for traveling step function. 


where I1(Z ,t ) is the finite-slope function shown in Fig. 27. 
Since 


2TH \= - Tish (L-a2'¢ 22 LU) casa) 
dz 

= res (o£2 ‘2 U- -02') (9 8b) 
then 

( 


L oot ul t 
a es rani, S66 =e 
J 3B (teyults %))az'= == dz = -Tlo) 


(99) 


Since for a step function I1(0) is also the current at the origin 


Ul 
0, then as before, 1(0) = 1(Z 20,t') = 


= 


I(t-r/c, and Eqs. 96, 
99) yield again Fq. (94). 


An alternative derivation of Eq. (94), using the for- 
mal mathematical properties of step and delta functions, proceeds 
directly from Eq. (81). 


A traveling step function, 
(95), 


using Eq. 
becomes 


T (2,t-Ric\= 0s 2. (\s6)) (100) 


The derivative in Eq. (81) ‘s 


—7 


T, ult- + - 2(-wse)) = T, J t- 0 - 2 U-esed) 


(101) 
where 5 is the standard delta function 


Eq. (81) becomes 


Hy = —-T _SM6_ sc 


ee eS ce Ae 
ra eee ae (t-wse@))d2 (197) 


Substitute* 


X= — eikrese) (103) 


at Cc: 
Aza- serena (104) 
“(58 


Eq. (102) becomes 


Hp = ttese Sa | SOCEM cos 
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‘3 L (ese) ¢ L-+¢ a) 
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- Oo othe (wise (106b) 


t 
The interval — 5 (iw 6) im aCe © can be rewritten as 


te + & f+ = (- (050) - Using h for L, Eq. (106) becomes 


Hat Axes _ WV(t-vley (107 


Sine 4irr 
Vie te v h 
—E TS wt Wes (107b) 
Cc Ca C 8) 
Thus the initial pulse of the radiation waveform ¢ Le step 


function excitation has the amplitude given by Eq. (106a) and 
lasts over a duration time interval of 4 (\-we) seconds, not 


seconds. This is discussed further below. 


“An alternative evaluation of the intepral in Eq. (102) is 
piven in Appendix VI. 


Since any function can be represented as a sum (super- 
position) of step functions, and Eq. (107) is the response of 
the antenna to a step function input, then by the linear super- 
position theorem (Ref. (33), p. i113), tne antenna response to a 


general current function input I, will be 


Hy lge) = Tex) rus aes 
4ur Sind 
Lene £4 (ro 1 
fete f+ SO) (2.08b) 


Eq. (108) is the same as Eq. (94). 

The radiated waveform for an input step functicen is 
thus given by Ea. (197). Actually, according *o the accelerated- 
charge radiation concerts, the sourzse of the radiation is in the 
moving leading edge (step) of this waveform. The distance R 
between this moving step and P varies, but for P sufficientiv 
far away, R may be replaced in the denominator by r in "a. (77), 
which is why r appears in the denominator of Ea. (107). 

The radiated field as observec at P commences at time 
r/c after the instant the step function enters the antenna at 
the oripin 0, and remains constant at the value given by £4. 
(1907) until the traveling current step wavefront reaches the an- 
tenna «end 7 * hy, where a reflected wave ‘s Initiatee, The Steld 


at P rema‘ns constant at the value viven ov fa, CL57) unet i the 


radiation from the reflected wave reaches :. Therefore the du- 
ration of this initial constant field at P is D (1-96) and not 


h/c. This can be seen dy examination of the following table and 


Fig. 28 
Observed time of 
Location of Local time arrival of radiation 

Event step wavefront at antenna at P 

1 entry at origin 0) r/c 

2 z Ore vic + ECi-ws6)\ (C 

3 reflection at h h/e ric + hCi-@se\ /C 

y return to origin 2h/e rile + 2hic 


In the local time frame at the antenna the timespan between 
events 3 and 1 is h/c, but to the observer at P, this time span 
is h/c (l-wse) » aS measured by his observations of the received 


s illustrated by the upper waveform in Fig. 29. 


he 


field. This 
The fact that the radiation-producinr sts. is travel- 
ing towards the observer produces the chservec time compression 
factor of (\- ws6e) . Similarly, the tire span between events & 
and 3 is again h/c in the local antenna tire, but ist h/c (le 090) 
as measured by the observer at P, This is because the radlating 


wavefront is now traveling away from the o»vserver, resulting in 


e . f om . 
a time exnarsson factor of C\+Ge9@). The time sran between 
events 4 ane 1 is ?h/e to both an onpserver 4° the antenna ond at 
P, as it should -°, cince both events cerur atio Fixed prolrt, 


the origin. 


The reflection of an incideni current wave I at the 
open end Z2 = h is equivalent to starting two waves from that 
end, each -I(t-h/c). One continues in the +Z direction (sup- 
posed to be extended) along with the incident wave, thus can- 
celing the incident wave for all Z2>h, and the other propagates 
back in the -Z direction (Ref. 29, p. 298), also (Ref. 30, p. 
105). The fields due to these two new waves are as follows 


(Ref, 30, pp. 104, 105 and Ref. 31, p. 144) 


+Z wove + Ny = =Tlt-We - vic) those" (109) 


47! sine’ 
a ates Hy —~-Tlt-hic- (A \-ese' (110) 
_ 4 r' sine! 


In Eqs. (109,110), nr! is the magnitude of the vector r’ drawn 
from the end Z = h to P, and ey is the angl+ between the +Z axis 
andr’. For P sufficiently far, r/ may be replaced by r in the 
denominators, and 6 by © . The radiation from these two waves 
thus appears to emanate from the end Z = h. 

To calculate the He field for step-function excita- 
tion of a center-fed dipole, or a monopole over a ground plane, 
it is necessary to take into account the radiation from both 
arms of the dipole, or the image in a ground plane. The vable 


below gives the observed magnitudes and time durations of the 


received radiation pulses due to both up and down-going current 
waves on each arm, over one time cycle of wave travel, starting 
with the step at the origin, followed by travel along the an- 
tenna, then perfect refiection at the ends Z2 = +h, and return 
to the origin where perfect absorption (no reflection) is assum- 
ed for simplicity. Fig. 29(a) shows the radiation wave-forms 
for Hy - multiplied by the common factor (4 v/To), The Z 
symbols in Fig. 29(a) show the directions of travel of the step 


on the antenna. 


Number in Direction of 4vr Duration of 
Arm Fig. 29 Wave travel Ha | =~) radiation pulse 


+Z AXXO get MM (i-036) 

Upper 1 (initial wave) Sin c 
~2 re WLesdete 2h 

Upper 2 (reflected wave) Sind c 

-Z A266, LNs 
Lower 3 (initial wave) Sino oct as 56) 

+2 ~ tore Nitwseyete 2h 
Lower 4 (reflected wave) sine val ce 


The radiation response due to an impulsive excitation 
may be found by taking the time derivative of the response due 
to a step excitation. Thus the radiation due to an impulse, as 
shown in Fig. 29(b), is obtained by time differentiation of the 
radiation wave-form shown in Tig. 29(a) due to a step. For@=90 


the two inner impulses in Fig. 29(b) ‘1erge at h/c to forma 


negative impulse with amplitude twice that of the two outer posi- 
tive impulses, as shown in Fig. 30(b). This impulsive response 
at 6=90, as well as for other values of @, can also be obtained 
from the response to a rectangular pulse of width T , by letting 
T approach zero. It should be kept in mind that these results 
assume no reflection at the origin. The case when there is re- 
flection at the origin is considered below. 

For B=50, the radiation waveforms of Fig. 29 reduce 
to those shown in Fig. 30. 

If the feeder transmission line impedance is not match- 
ed to ihe surge impedance of the antenna, then the step wave- 
tronts are reflected on return to the origin Z = 0. Let the 
current reflection coefficient at the origin be Ko (° 0.72 in 
Section 3’. Further, for the sake of completeness, assume non- 
perfect reflection at the ends Z = +h, and let the reflection 
coefficients there be Ke (= -0.9 in Section 3). Then the radia- 
tion waveforms shown in Fig. 29 must be modified in amplitude 
and extended in time as shown in Fig. 31. For the sake of defi- 


niteness, @= 45 ,o = 0.72, and We = -0.9 in Fig. 31. 


4.3 RADIATION DUE TO PULSE EXCITATION; COMPARISON WITH A 
FREQUENCY-DOMAIN ANALYSIS 
In this subsection, the radiation due to a rectangu- 


lar pulse will be calculated using the time domain techniques 
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FIG. 31 Radiation H,(atr/T.) due to step function on a dipole. 


discussed in the previous subsections. Let the pulse duration 
be YT seconds, where \ is less than the travel time h/c on either 
arm of length h of a center-fed dipole as shown in Fig. 24, or 
on the length h of a monopole over a ground plane. Such a pulse 
traveling on each arm of the antenna may be represented before 
reflection at the ends as the sum of two traveling step func- 


tions as follows: 


U (t- Zic)- U-T- 2ic) (llia) 
UY (etzle) -U t-Tt2l¢) (111d) 


upper arm: 
lower arm: 
The time history of the radiation waveforms due to the 

first terms (leading edge of the pulse) in Eqs. (llla,l11b), in- 
cluding reflections at the ends Z = + h and at the origin Z = 0 
has been illustrated in Fig. 31. The radiation waveform due ‘> 
the second term (trailing edge of the pulse) has a similar shape, 
except that the sign is revised due to multiplication by -1, and 
the waveform is shifted to the right along the time axis by T 4 
due to the time delay of JT seconds between leading and trailing 
edges. Addition of the two sets of waveforms thus gives the 
radiation waveforms for rectangular puise excitation of a 
center-fed dipole, with a current reflection coefficient Ke at 
the dipole ends, and we at the input (2 = 0), for a pulse width 


less than the antenna travel time h/c. 


The details of this time-domain calculation have been 
carried out, but are not repeated here, to permit cor -irison 
with an example given in Ref. (34). In Ref. (34), the center- 
driven cylindrical antenna has a length to radius ratio he 504, 
with a source excitation ulse of trapezoidal shape with a base 
width of about 0.55 h/c, a flat top width of about 0.26 h/c, 
and a width ot about 6.38 h/c at the 50 percent peak amplitude 
point. The transmission line has a 50 ohm impedance. To defin- 
itize the comparison, the following numerical values were assum- 


ed: 
vy20.4hle 


The patterns in Ref. (34) are calculi ted for the far- 
zone field using an inverse Fourier transform of the response to 
jut 


sinusoidal excitation C . A comparison between these patterns 


as given in Ref. (34), and those calculated as explained above, 


assuming a rectangular pulse of wiuth Y = 0.u h/ec are presented 


in Fig. 32, for different values of polar angle 8 . The patterns 
of Ref. (34), as shown in Fig. 3° nave been simplified ‘u that 

's small amplitude oscillations nave been omittea, however the indi- 
cated amplitudes, widths and tires of the princisal pulses follow 


the values piven in Fef. (38) 2 closely as can be read trom tne 
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curves presented therein. The time-domain results have been 
Sc 

normalized so that the first pulse for 6290 (broadside) has 
the same amplitude (=1.£5) as the corresponding pulse in Ref. 
(34), thereby afforaing ready comparison. 

Inspection of the two sets of curves in Fig. 32 shows 
a fairly reasonable agreement, considering the relatively ap- 
proximate nature of the ime-domain calculations, except as 9 


e ° o ° 


approacnes small angles. For @ = 30 , 70 , 60 , and 45 , the 
agreement is fairly rood as measured by the number of pulses, 
their amplitudes, widths, and positions. For example, as & 
decreases, the amplitude of the first pulse increases approxi- 
Mately as predicted by the (14c036)/ Sine factor of Eq. (107a), 
until 8-20. Further, using time-domain analysis, it is simple 
to show tnat us @ decreases from 30 , the received pulse width 
remains constant at Tountil i.e anple @ is reached at which the 
quantity (\-Cos@)n/c becomes less than T. For O less than this 
angie, the pulse width, tor an ideal rectangular pulse, becomes 


C\-cas@)ii/c. lor VY = u.4 h/eo, this anple is 53. Inspection of 


the curves piven in Ref. (34) shows that the received pulse 


° 

widcu bLewans to decrease for ® around $5 , in apreement with the 
tire-domsin prediction, and continues to decrease as @ cecreuses, 
agyate in gyreertent with the tive-dorain analysis. 
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However, it is clear that there is serious disagree- 
° e 
ment between the two sets of results for 9@€26 . The time- 
domain analysis as presented in this report, which leads to eq- 
uations such as Eqs. (70) and (107), which agree with previously 
published work (Refs. 29-32 incl.), predicts a field which con- 
tinues to increase with decreasing 6 as 


lim _ \4eosq ae ce (1.13) 
G--Ssmall sine (>) 


On the other hand, while the frequency-domain analysis of Ref. 
(34), presents curves which show a field which does increase 
approximately as (\tcog6 ) /sine, from 90° to "0 » the B=20 
curve distinctly shows a field which has begun to decrease again. 
The disagreement arises in the singular behaviour predicted by 
Eqs. (70,107,113) for the field which increases without limit as 
© approaches zero. This behaviour is discussed in more detail 
in subsection 4.4 below. 

However, it is stated in advance that while this dis- 
cussion shows that the total radiated energy of the pulsed di- 
pole antenna remains finite, and that the fields do not go to 
infinity as Q~+O for an infinite-length antenna, nevertheless 
the disagreement noted above for small values of © remains 
unresolved. pxrperimental investiyation of the fields of a 


pulsed dipole for 8 = small might help to resolve this point. 


4.4 DISCUSSION OF PULSED WIRE-ANTENNA RADIATION, ESPECIALLY AT 
SMALL VALUES OF 6 

Eqs. (70,71) have been derived in this report for a 
current wave traveling in the +2 direction on a finite-length 
linear antenna, such as a dipole, the effect of the ends not be- 
ing included. The same equations are derived using the retarded 
vector potential A for a finite-length antenna by Ross et al 
(Ref. 31, pp. 17, 117) and also for an infinite-length antenna 
using a Green's function by Ross et al (Ref. 31, p. 142). The 
effect of the ends is accounted for by a reflected current wave 
by Ross et al (Ref. 32, p. 115). The same equations were pre- 
sented without derivation by Manneback (Ref. 29, p. 294) for an 
infinite-length wire; he showed that these expressions satisfy 
Maxwell's equations, neglecting \/ r* terms compared rol /¥ terms, 
and that for © smail, Eqs. (70,71) reduce to the expected fields 
as We of a plane electromagnetic wave traveling at velocity c 
along a wire. Manneback also took into account end effects by 
partial or total reflection at the ends; see Eqs. (109,110). 
Schelkunoff (Ref. 30) shows that Eqs. (70,71) are really exact, 


2 
including \/ Y terms, for an infinite-length wire; he also 


accounts for end effects by equations such as Eqs. (109,110). 
es The purpose of the above 4iscussion is to show tiat 


ly * the radiation equations for finite-lenpth antennas are the same 


as for the infinite-length antennas, except that provision must 
be made for end effects. For infinite-length antennas there are 
no end discontinuities; the fields in the direction of the an- 
tenna (small@) travel with the current wave along the antenna 
at velocity c. For the finite-length case, the ends represent 
a discontinuity, such as partial or total reflection, or absorp- 
tion. The radiation from the finite-length antenna is thus cal- 
culated using the same equations as for the infinite-length 
case, modified by appropriate reflection coefficients for termin- 
al conditions, with due regard to the direction of current wave 
travel and time delays (e.g., Eqs. (109,110)). 

For an infinite-length antenna, Eqs. (70,71) present 
no special problems as 87%O. Ross et al (Ref. 31, p. 137) show 


that Ampere's law in the form 


p eye ee L (114) 


is satisfied by Eq. (70) for 870. See Fig. 33. Thus 


qa y sin6 
& HED) -32 a Howard = H cet) 2 itr sine 


— Llt-rle) tess | anvane 
47 ane 


Using \+¥@$6 7 2 and Y-e2 as OO , 


® Ge, ty. dl = U(t- 2lc) (115) 


Thus the line integral of H near the wire at distance Z from the 
origin, at time t, is the current that was at the origin at time 
Z/c ago. Since the current is assumed to travel undispersed 
along t'e wire at velocity c, this is exactly the current that 
is now (at time t) at 2, and has reached the loop area. It is 
pointed out thac Eq. (114) is not the complete form of Ampere's 
law, which is, for stationary boundaries, 


g Fal = T-85 + & Se] E 45 (116) 
5 


The first term in the RHS of Eq. (116) is the conduction current, 
and the second term is the displacement current. However, as 
Q@--Q, and the surface of the wire is approached, 33 + AS as - 
the radial component (Z-compone t) of £-»9, and there foref-d520. 
Thus Eqs. (114,115) hold as O70. Schelkunoff (Ref 30, p. 109) 
arrives at the same results using a conductor which is a semi- 
infinite cone of half-angle W . As 8+) (surface vf the cone), 
Ampere's law in the form of Eq (114) is invoked, the radial com- 


ponent of E being zero at the surface. ‘See Fiy. 3b, 


FIG. 33 Line integral of H for semi-infinite wire. 


ric. 34 Line integral of H for semi-intinite cone. 
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The sore2going discussion indicates that as far as cal- 
culation of the fields due to traveling waves on infinite wires 
is concerned, including the case of 96-0, there is agreement 
among the various writers cited. Fig. 35 illustrates how the 
fields wavefront travels with the current wave along the antenna 
at velocity c. Two wavefronts are shown at time t. The lead- 
ing wivefront is due to a current element I, 42 now at Z= NY, 
and formerly at the origin at time t-Nic , the second wavefront 
is due to a current element eae now at 2=X2 , and formerly 
at the origin at time t-telc. 


For a finite-length antenni, Fiy. 36 illustrates how 


the radiation field at a far-field ;oint Ff varies in time. As- 
sume a step-function excitation. ihe initial radiation to reach 
P is due to the current wave traveling in tne +7 directicn. At 


tine i (i @s6) 1 ¢ later, the radiation from the reflected cur- 
rent wave arrives at P. At exactly 8=O, the two wavefronts are 
just tangent. Calculation of the resultant field at Q@20 is 
beset with a number of difficulties such as the boundary condi- 
tions at the end © = h, the effects of a finite wire diameter, 


the discontinuity experienced by the traveling field wavefront 


as it moves off the end of the wire into Space, and the zero 
duration of tnis field. however, at non-cere valies o1 @, it 
is clear that the radiation tnat arrives first at } ts due only 


| vadiation wavefront 


vadia tion wave- 


front due to T, 


FIG. 35 Currents and wavefronts for current wave 
on semi-infinite wire. 
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to the initial current wave traveling in the +Z direction; the 
radiation change due to arrival of the current wave at the endh 
reaches P at time h (i-cos@)/¢ later. Hence during the time 
interval Oft<hl-os6)/c, the curreni step has not yet reached 
the antenna end, and the radiation at P must be independent of 
what occurs at Z = h. Therefore, during this non-zero time in- 
terval, unaffected as yet by end conditions, the field is given 
by Eq. (107a) for a step waveform, and by Eq. (108a) for a gen- 
eral current wavetorm. All the cited references (Ref. 29-32 
incl.) agree on this form of the radiation equation. Hence the 
bothersome question remains as to the behaviour of this field 
for small © as 2 [sine » which disagrees with the calculations 
of Ref. (34). 

In terms of field amplitude or power per unit area, 
the linear pulsed antenna appears to act as an end-fire antenna. 
Ordinary (sinusoidal time variation) traveling wave linear 
antennas also act to tilt the radiation pattern in the forward 
direction, the tilt increasinr with the ratio of antetna leneth 
to wavelength. However, such antennas are terminated at the 
ends so as to eliminate standing waves, i.e., an absorption 
termination instead of an open end. Hence any comparison with 


ordinary traveling wave antennas must be t’: ated with cauticn. 


In terms of radiated energy, some interesting conclu- 
sions are easily arrived at for tie pulsed dipole antenna. As- 
sume step excitation, for simplicity, and perfect end reflec- 
tion. For one cycle of current travel on the antenna Oft£ 2hic, 
Fig. 29(a) shows the radiated field, for 627 /2. ror See ; 
it is simple to show that the total radiation waveform remains 
unchanged in form, consisting of two pulses as shown of tne same 
amplitude and time duration as for G@lW2. nence the radiated 


energy per unit area (two pulses) at any angle Qis 
ahic 5 ae 
E(®)=| jute Wdta=e2 jt —4o___h (\-\esal) 
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& (4m sme < 
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o\—| A to | (117) 
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E (0) ~ (118) 
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Thus the energy per unit area is tinite at all anyles. Also the 
F - ° : ; : 

broadside (= 90°) energy per unit areu is twice that in the end- 
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fire directions (Q@=0 , \BO decreasinye Smootily from tae 
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and power as 9-~ 0,180 is not manifested in the energy. In 
terms of radiated energy, the antenna does not act as an end- 
fire type. 

The total energy cx radiated in cone cycle on tszhie 
by a step excitation on a dipole antenna is found by integrating 


Eq. (117) over the surface of a sphere. This leads to 


2 en we 
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Te = + 2 Ww2 
C2. 2G (119) 


Hence the totai radiated energy is finite. 

Before leaving the subject of the singular behaviour 
of the field amplitudes as @-~O, it is pertinent to note what 
happens if the velocity of propagation Vo of the traveling cur- 
rent wave along the antenna is very slightly less than c. Then 
the factor of C\+c03 8) (Sin = 31n8 (Cl-aae) in Eq. (107a) is re- 
placed by psme/Ci-gcose), where GEUp/C<} (Rei. 32, p. 123). 


The significance of this is that the "infinity catastrophe" at 
6-0 disappears. The field now goes to zero at @=O , and has 

a maximum at OQ, =¢05' 8 given by P/fsmO, . For example, for 
$20.98, Om =.S° . It is only a conjecture if some ef- 
fects of this nature are indeed taking place on bare wire an- 
tennas. The discussion of propagation of sharp wavefronts on 
lossy lines in Appendix III indicates that the definition of an 
effective velocity of propagation depends upon the details of 


build-up of the wavefront at a point along the antenna. 


SECTION 5 


RADIATION FROM LOOP ANTENNAS 


5.1 INTRODUCTION 

In this section, the radiation of a circular loop an- 
tenna ‘s examined. Two examples are given. The first is a loop 
driven with sinusoidal time-varying excitation, for which the 
solucion is well-known. The radiation field is derived from the 
accelerated charge Eqs. (V-7,32), which for a linear element dZ 
lead to Eq. (74). This example illustrates that these equations 
produce correct results for known cases. The second example is 
a loop driven by pulsed excitation, for which a solution has re- 
cently been made available using both frequency and time-domain 
analysis (Ref. 35). The time-domain analysis is examined to 
verify that its formulation in the time domain agrees with that 


developed in this report. 


5.2 SINUSOIDAL TIME-VARYING EXCITATION 

The loop is taken in the XY-plane, with radius b suffi- 
ciently small compared to wavelength so that the current can be 
assumed to be constant around the loop, to a good approximation 
(Ref. 22, p. 56). In Appendix VY it is shown that Eq (V-7), for 
the radiation field of an accelerated charge moving at low relat- 


ivistic velocities, becomes Eq. (V-32). For a loop, differential 


oe ee 


length dZ is replaced by the differential circumfe* ential length 
dl dL ay = bdo’ ay’ (120) 

See Fig. 37. The coordinates of the current scuce element dL 

are ( b,w/2 4! ), and those of the field pé*.t P are (1, 8,9). 


Eq. (V-32) becomes 


AW Tdoixt (121) 
41¢ IR 
vet 


For point P in the far field, R may be replé.ed by yr in the de- 


nominator of Eq. (121). Using 


7 ited 2 T (t- Ris (122) 
ret 
T=TI, oist (123) 


it follows that 


IT (t-RIC)= pore ™e i 
at 


(124) 


Hence 


as = ihe ae ‘ 
au — wt.’ = Sux h eset (125) 


For P in the far field, R is essentially parallel to Yr, so unit 


A n 
vector Y\ may be replaced by unit vector Gr , and 


R » V- b ws (126) 
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FIG, 37 Loop geometry. 
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where, from spherical trigonometry 


CoS) = WS (b-¢°) 30 (127) 


Hence 


“KER SAKE jebcos)  —jer jeboo(d-6') sine (128) 
ce =C og =e e 


by hypothesis, «tb = aub/dheel | hence 
@ ikR ~ @ XT (it jb cos c-6"\sine (123) 
Resolving dL along the directions ay parallel to Pf, the pro- 
jection of r on the XY-plane, and ay (=a4) perpendicular to ©, 
diadt [‘os (o-¢') ay rstn (o-6') a, | (130) 


Substituting Eqs. (120,129,130) into Eq. (125) and integrating 


: ‘ 
over the variable o ; 


ew 
H = {wt © IS EIA, | [ieyebamne cod (a-¢')| 


4ycr 
foa 
° [cos (4-4) a ;* sm (b-4)4, |b Ka, 


4 --\wt. b EWC (tebe) 4 Lx4y) 
4ner 


2 
Using loop area A= Th ’ r= w/c , Magnetic moment Mi, = wa > 


-~ ey - 


. 3 (wt- Kr) be 


¢ 
: 2 
G= —TMelk™ sino (131) 


which is a well-known result (Ref. 21, p. 95). 


5.3 PULSE EXCITATION 

In Ref. (35) the loop is as shown in Fig. 37. The E 
field is calculated in (Ref. 35) using the radiation field of an 
infinitesimal dipole given by Sommerfeld. Using the equations 
developed in this report, E is given by Eqs. (II-34, V-37) as 


f= fare (Uxi) (132) 


Hence Eq. (121) can be used, to yield 


ab -_ \ (vated) xi (133) 
4atec® R 
wh = es 
From Eq. (V-22) for the differential current element hn xdl , 
ya ° _ 
2) BA-Ric) = Tt-Ricyab (138) 
at’ 


In the far field R is replaced by r (except in the arguments of 
Dp 7; Dp y N - =m A 
P or IT), and since R and r are parallel, then N=AV/V=S Gy , 


Hence Eq. (133) can be written as 


2 = ~ 
a =—t_ SB Pet red xF bx ca) 
antec’ es } at® 


Eq. (135) is identical to the equation for dE developed in (Ref. 


35, p. 28). 


aren OT] 
n 
. 


ee Ata, 


= yee 


In Ref. (35) the source and field point coordinates are(b, £, o,\ 


and CR , do) respectively, whereas here they are (b, Ez, b') 


and Cv) 6, d) . The distance bstween source and field points is 
r in Ref. (35), whereas here it is R. Taking this into account, 
Eq. (135) is the same as that used in Ref. (35) (p. 28, Eq. (4)). 
rN A 
Returning to Eq. (133), using Eq. (130), AXA 24 F 
A A 
WW Qa, , there is obtained 


4 (aL xf )af =-ab |coatd-¢') Ay tsinld-d 96 ae | 


(136) 


From Eq. (123), dE has two components 
ag. --___|. ss A Tlt- Ric) dLsm o-4") 098 1 374) 
eo  4Tec®r at 


aby =. \ BD Tt-Ric\ dl. os (4-4") (137b) 
4iec@r at 


In Ref. (35), the excitation is a pulse applied to the 
input terminals shown in Fig. 37. Similar to the case of the 
linear dipole discussed in subsection 4.3, this produces a cur- 
rent pulse traveling circumferentially on each arm of the loop, 
in opposite d directions. Each pulse is represented, as in 


subsection 4.3, as the sum of two separated step functions, re- 


sulting in four step functions in all, It will suffice, for the 
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purposes of this report, to investigate the radiation due te one 
step function; the detailed results of the superposition of the 
four step functions is given in Ref. 35. 
A step function of current traveling at linear veloc- 
ity c on a loop is given by an equation similar to Eq. (95), 
( 
with Z replaced by circumferential length L= bh - Thus 


DT. UCe- be) (138) 


Using Eqs. (126,127,133) 


2 — y- bcos (d-6')aind (139) 
Tt-R()=T,U(t- F-2[o> wage }sne|) caso) 


1, OTH RA) = T, §fe- £-2 [b-wa1d-4") 108) 
ot (142) 


( 
Substituting Eq. (141) into Eq. (137), using d= badd » and 


t 
integrating over the variable d gives 


ev ( 
ae bcos 6 To (b-¢' $.L4€(o'}) ae (142a) 
= bane (Goce dtd : 


(o) 


2 ‘ ‘ 
jog Be[ Leet 


(2) 
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F(¢\S- 24-5 (d-¢)ome] G4) 


The integrals in Eq. (142) are of the type discussed in Appendix 


vi, 
Jaa \slFev-a}de = pated (au) 
ad! | 


t 
Here a; is a zero of the argument of the delta function, i.e., 


: :. Lie v=d . From Eq. (243), 


A £14) = & \- sin (6-4! \sine (145) 
‘ ad’ d: C 


Hence Eqs. (142) become 


_. -heto28 Sin (o-$% \ (146 a) 
4Tlecr Ss ain (o-4,'\ sine | 


' Eg= _ Le tortd-d.>) { (1"6b) 


4Tecr | I-om (o-4/\sine | 


Tl 
® 
i 


Eqs. (146) are the same as those in Ref. (35) (p. 31,32, Eqs. 
(12,13)), using vi (traveling-wave velocity) = c, except that 
there is a + sign in front of the sin (O-4; \ sme term. This 


may be traced back to what appears to be a typographical or a 


sign error in Ref. (35). If one substitutes the equation 


¥= R, —b cos Cb-d,) Sin 6, (Ref. 35, p. 30) 


into the argument of the delta function 


bdo 
(4-T, -28-+t (Ref. 35, p. 30,31 
J . We oD Eqs. (10,11)) 


one obtains the argument 


+ 7; - be 2 +s te £005 (ded) ne, 


whereas a — sign appears in front of the Co (4-42) in, term 
in. Ref, ®5 {p. 31). 
2. is concluded that the time-domain analysis of the 


puized loop problem in Ref. (35) and the analysis developed in 


this. report are very similar to each other in basic formulation. 


SRCTION 6 


TRANSIENT RADIATION FROM APERTURE ANTENNAS 


1, INTRODUCTION 

The radiation from aperture antennas such as horas and 
paraboloids is calculated in many analyses from diffraction of 
the aperture fields Es é Hh. For reflector types the reflec- 
tor-surface fields or currents can also be used (Ref. 21, p. 
144). For wire antennas the preceding sections have shown that 
the use of accelerated-charge radiation leads to a time-domain 
formulation in terms of the time-derivative of the currents on 
the antennas. A similar analysis for reflector-aperture an- 
tennas, using the accelerated charges on the actual conducting 
surfaces to cal..ulate the transient radiation has not been at- 
tempted as yet. Instead, in this report, a general time-domain 
method of calculating the transient radiation from aperture an- 
tennas developed by Chernousov (Ref. 37) is reviewed, and some 
important results are neted. This method uses the antenna sur- 
face or the aperture fields and the Huygens-Kirchhoff principle 
to calculate the radiated fields. It is shown in this report 
that the radiation from the equivalent current sheets which re- 


place the aperture fields in Chernousov's analysis, may also be 


obtained from equivalent acceleracing charge sheets over the 


aperture. It is also shown that Chernousov's results for the 
radiation of a planar in-phase aperture in principal planes 
reduce to the results for a one-dimensional antenna given by 
Cheng and Tseng (Ref. 38). Some recent work by Maddox ‘Ref. 3) 


is also described. 


2. REVIEW OF CHERNOUSOV PAPER ( JF. 37) 
Consider the antenna aperture fields E (t,t), WCts yt) ; 
which must be known, and their equivalent currents, which are 


arbitrary functions of time, 


Ke (%,¢) = mM x KCB E) (147) 
Ky Cs, 4) =-W XK EC t) (14 7b) 


where Ce Su are the usual equivalent ficticious electric-and 
magnetic-current sheet densities over aperture S (Ref. 6, p. 


486). The retarded potentials are (Ref. 6, Chap. 13) 


A, Ct ey oe “| Ke(G,A RIV) Js (148a) 


Ayce a) =e 


Ku Ce k- 80) CX. t- Rls) AS (148b) 
41), R 


See Fig. 38. Ys, is the vector from origin 0 to aperture area 


element a3 » and v is the propagation velocity in the medium. 


a a 


Lo 


aperture S 


FIG. 38 Apert re antenna. 
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Then, following well-known methods, the fields Echa, H¢7,t) 
at P are expressed in tarms of vector operations on Re and Ba . 
The results, which are not repeated here, are given in Ref. (37), 
and hold for near-as well as far-fields. 

Chernousov defines the far field ty the following 


equations, 


— 


1 |_2&¢, Be ale ae : 
eo le >? | (8,6) (149a) 


Alea, ‘)| >? 5 (Fee) at | (149b) 


plus the requirement that R22? linear dimensions of the antenna, 
where F denotes E or H. To illustrate these equations, for con- 
ventional sinusoidal time-varyirz excitation, let Fo €(H jest . 
Using WuC , Eq. (149) reduces to Re Nazi » Which is a well- 
known criterion for the far field for eset excitation. For 
arbitrary time variation, using Chernousov's results (Ref. 37, 


Eq. (8)) and Eq. (147), the radiation field can be shown to be 


- = -_ ae t Teeth 
WCF = ae & Re (i,t + et) 


47cr 
5 a p 
A te. fe. aly Nes 


jaye 


(150) 
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—— 


‘ A 
where C =t- yl, ana ¥Y is a unit vector along Y . 


3. ACCELERATED-CHARGE VIFWPOINT 

It will now be shown that Eq. (150) is consistent with 
the accelerated-charge radiation results developed in Eqs. (68, 
74, V-7). For a one-dimensional current flow along the Z-axis, 


these equations show that 


Ah HS pas pene dS T(4t- E)\dzxR (152) 
4T1CcR dt 


a 
where R is a unit vector along R, the vector from dZ to field 


point P. In a similar fashion to Eq. (68), which replaces the 
charge-acceleration product by a product of the time derivative 
of the current and az » for a surface sheet current flow Ke 


amps/aeter, the charge-acceleration product would be replaced 


by the following product: 


Lae] = — Ke Cran t- &\ds (152) 


The: the radiation due to the equivalent (ficticious) electric 

current sheet Wig which replaces the true aperture field WC ED 
can equally well be considered as arising from the radiation of 
equivalent (ficticious) accelerated charges as given by Eq. 


(152), leading to the following two-dimensional generalization 


of Eq. (151): 


- SW - 


ab (¥t) - —1__ Deh tIe)SSK RK case) 


47 cR ot 
For the far field, RX - Yet , R is parallel to r, and Eq. (153) 
becomes 
a A = = r Yoo 5 
~*~” 1 %.-4 
AW (Ft) a- Cin Bh GG t- f+ Ft) sy) 


4c ot 


Eq. (154) is identical to the Ke term in Eq. (150) from Chernou- 
sov's paper. A corresponding term arises from the Ky current; 
this is the Ky term in Eq. (150). Thus the Chernousov results 
for aperture fields can be interpreted in terms of equivalent 
accelerating charges. 

In the case of linear antennas, the charge-accelera-~ 
tion product was formulated in terms of the time derivative of 
the actual current. For an aperture antenna, the above shows 
that a similar correspondence obtains between the time deriva- 
tives of the equivalent (ficticious) aperture currents and ficti- 
aperture charge acceleration products. A transient radiation 
analysis in terms of the actual reflector currents has not been 
attempted, but presumably would produce a similar correspondence 
between these current time-derivatives and the actual surface 


charge-acceleration products. 


4. PLANAR IN-PHASE APERTURE 

Chern>rusov examines transient radiation from a planar, 
rectangular aperture excited in-phase such that the aperture 
field at any time has the same value at all points in the aper- 
ture. Thus the aperture excitation may be written as 

Es(¥, +) = E,lé €,C%) (155) 

where e, (Te) may be equated to unity over the aperture. The 
aperture is taken in the XY-plane, radiating in the Z-direction, 
of dimensions ~@/(24¥ © al2,-b/a¢ ae b(2. The aperture field 


— 


E 


5 is polarized in the Y direction, and BS is in the -X direc- 
tion. For the case where Es, Ws form the front of a free-space 
TEM wave, Chernousov shows that the radiation pattern in a prin- 


cipal plane such as XZ is of the form 


EG) = Ab Litese) ok fest ¥2\- Be ‘a 
mee oo (156) 
1 


where Q@ is the usual polar angle with respect to the Z-axis, 
t'= t-tle , and Ta = asme@/(2c). Since Vy is the time dif- 
ference between radiation from the aperture edges and the origin 
it is seen that the radiation is the sum of two waves, which 
appear to originate at the aperture ends, and have the same time 
form as the excitation field. 

Cheng & Tseng (Ref. 38) consider a one-dimensional an- 


tenna with current eyvcitation C(¥,t\ » When this function is 
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separable then UY, e\ = Aly) f(t) - For a comparable example 
with the above Chernousov case, select the example in Ref. (38) 
where AWN=\) for -al2¢xSal2. For @sQ ¥I (Ref. 38, 

footnote 3) so that (letose\ ¥ 2 , the radiated field is (Ref. 


38, Eq. (13)) 


E(6,4\ = de _}f Ctl tu)-4({e]-at)} (157) 


( 
where You asm /C2c\= Te » and {t]=t- Vicxt - Hence 
apart from differing amplitude factors, which is expected since 
one analysis is based on aperture fields and the other on cur- 


rents, the two patterns of Eqs. (156,157) are identical. 


5. SOME ADDITIONAL RESULTS 

Using Chernousov's formulation, Maddox (Ref. 3) has 
investigated the scattered far-field due to a conducting disk 
whose surface is perpendicular to an incident pulse-type plane 
wave. In the broadside direction (6 =o) » the time waveform of 
this field is shown to be the negative time-derivative of the 
incident waveform. Thus the response to an impulse illumination 
is an inverted doublet. This has been verified experimentally 
(Ref. 39). Maddox (Ref. 3) has also derived the equations for 
the radiation E field on the axis of a paraboloid antenna illu- 
minated by a pulse-type sperical wave whose orig'\n is at the 


focus, again using the Chernousov formulation. Calculations and 
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interpretation of the results are not yet complete; however it 
appears, depending upon the geometry of the problem, that it is 
not always true that the radiated field is the negative time 
derivative of the incident excitation. 

A limited amount of experimentation has been done, 
using available equipment. The transmitter is a Spencer-Kennedy 
Laboratory pulse generator (mercury switch), with a one nano- 
second pulse width (one foot resolution) and a rise time of 
about 350 picoseconds. This excites a long vertical monopole 
from a coaxial line. For a receiving antenna a linear string 
of five carbon resistors has been used (Ref. 26). The display 
oscilloscope is aH. P. 185B with a 187B preamplifier. With 
this equipment in a non-anechoic roon, measurements have been 
made of the returns from a parabola (3 ft. diameter, 9 in. focal 
length) and two disks (1, 2 ft. diameter). Qualitative agree- 
Ment has been obtained with published results (Ref. 40). More 
refined experiments on aperture and wire antennas would require 


test equipment with considerably higher resolution. 


SECTION 7 


CONCLUSIONS AND RECOMMENDATIONS 


Based upon the results presented in this final report, 
the problem of radiation from antennas excited by impulsive sig- 
nals can be approached in the time domain from the viewpoint of 
radiation from accelerated charges. The transition to antennas 


is accomplished through derivation of a radiation equatjen which 


expresses the accelerated-charge radiation in terms of the time- 
sz derivative of the antenna currents, either real or equivalent. 
Applied to pulse-excited dipoles, linear antennas and aperture 


antennas, this radiation equation produces a result derived pre- 


viously by retarded potential and other methods. For a pulse- 

excited loop, the result is the same as that obtained by another 
a investigator using the Sommerfeld radiation equation for an in- 
finitesimal dipole. For sinusoidal time excitation, the above 
radiation equation reproduces well-known results for standard 


antennas such as small dipoles and loops, a half-wave length 


a dipole, and a linear traveling-wave antenna. 


It has been shown that the radiation characteristics 


4 of antennas may be expressed in the time domain in a number of 


equivalent ways, involving the integration of either the retarded 


potentials, or the Sommerfeld dipole radiation formula, or the 


radiation from accelerated charges over the differential elements 
of the antenna. If one accepts the premise, as preferred by 

this writer, that the basic physical radiation mechanism is the 
presence of accelerated charges, then the other methods may be 
interpreted as alternative mathematical formulations. 

In this report, the radiation fields due to acceler- 
ating charges in all cases have been expressed in terms of the 
time-derivative of the actual or equivalent currents on the an- 
tenna, which are assumed to be known. Then the method of rétard- 
ed potentials, for example, could also be used, without any: ref- 
erence to accelerating charges. Thus, based upon the work com- 
pleted thus far, it is concluded that while the accelerated 
charge radiation approach provides, in the time domain, a direct 
physical explanation of and an analytic basis for impulsive an- 
tenna radiation, it has not been demonstrated that it has sig- 
nificant advantages compared to other methods which also use 
known or assumed antenna currents in their formulation. It is 
probable that this conclusion will be altered by the results of 
further work on transient antennas with new configurations, where 
the currents or aperture fields are unknown and must be solved 
for (i.e., time-domain boundary value problems), or must be con- 


trolled by new and novel techniques. 
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In the report, the radiation from a pulse-excited di- 
pole is examined at various angles, using a time-domain equation 
which has also been derived by others. Comparison with a pub- 
lished frequency-domain analysis shows reasonably good agreement 
at all angles except for angles close to the end-fire direction. 
Further work on this point is recommended. Experimental data on 
the radiation from a pulsed dipole (or monopole) at all angles is 
also recommended. 

The receiving characteristics of antennas in incident 
transient fields have not been investigated in this report. It 
is noted that the receiving response of a monopole over a ground 
plane, matched at the base, to an incident impulsive plane wave- 
front, as given by Ross (Ref. 27) has a waveform identical to 
that for the transmitted radiation of a step function given in 
Fig. 29(a) (Ross's 8 is the complement of the ® used in this 
report). Investigation of receiving characteristics and reci- 
procity relationships with transmitting characteristics is an 
obvious recommendation. 

A promising method due to Chernousov for analysis of 
the performance of transient high-gain aperture-type antennas 


has been identified and described. This is a gener: oach, 


based upon replacement of known or assumed apertv by 


equivslent surface current sheets. It is shown in this report 
that the Chernousov results can be recast alternatively in terms 
of equivalent surface accelerating charge sheets. For a one- 
dimensional case, it is shown the’ the Chernousov formulation 
reduces to results derived previously by Cheng and Tseng. It is 
recommended that the analysis of impulsive high-gain aperture- 
type antennas be extended both in scope and in detail, including 


primary feed antennas, to further the state-of-the-art of this 


important class of antennas. 


12. 
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APPENDIX I 


LIENARD-WIECHERT POTENTIALS 


1. SCALAR AND VECTOR POTENTIALS 

The Lienard-Wiechert potentials are the potentials result- 
ing from charges with arbitrary spatial and temporal distri- 
butions. From toese potentials the E and H fields due to 
accelerating charges may be found. For the sake of completen:- 
a derivation of these potentials is included in this Appendix. 


Maxwell's equations in Gaussian units (Ref. 1) are 


wxi= 403+5 28 os 


Wit =-2 28 (1-2) 


VW DO =4T7P (1-3) 

V8 =C) (I-4) 
The constitutive equations are 

D=e.e. B=Ar:s Js (1-5) 
Because of Eq. (1-4) and the vector identity div curl = 0, B 
can be derived from the curl of a vector potential A: 
B= Wr A (I-6a) 

Similarly, because curl grad = 0, E may be aerived from a 


scalar potential $ and A by 


E=-VG -1. OB (1-6b) 
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Thus the curl of Eq. (I-6b) recovers Eq. (I-2). The Lorentz 
condition (Ref. 1, p. 180) defines div A by 

yh--- 28 (1-7) 
Vector manipulation of the above equations leads to the 


standard wave equations 


w*6 - ae s-arp (1-8) 
q¢n - L 28. 407 (1-9) 


where J is the impressed current density due to external sources. 
For static fields Eqs. (I-8,9) reduce to the Poisson 
equations 
W'$=-4iP (1-10) 
27-407 (I-11) 
WA =- 409 


The solutions of the static Eqs. (I-10,11) are the particular 


integrals 


{ 
. 


OCk \ = \ oe az! (I-12) 
AGQ)= | Ay a dx! (I-13) 


Here x is the position vector of the field point P(x, , Xp» Xs ); 


— { ¢ 
x! is the position vector of the source volume ax! = AXidxy dX, 
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‘ ‘ 
at point P( X% ) Xe’, X35 ), and R is the magnitude of the 
vector R= X-x' from e toP. See Fig. 39 
For dynamic fields, the time-dependent Eqs. (I-8,9) are of 


the form 


TY- =p ay =-47 £(%t) (1-14) 


The time-retarded solution of Eq. (I-14) is (Ref. 1, p. 186) 


Yat =| Alt-ttRle) CL) a ‘dt! ca-1s) 


where S is the Dirac delta function. The integration over t' : 
using the properties of the delta function 

§ (ti -t+Ric)=0 oe 
except at 


t'=t- Ric (I-17) 


mtg 
Yas | L8G) tleet ay (1-18) 


yields 


t 
where [ ace means that { is the "retarded time" given by 
Eq. (I-17). The time-varying solution Eq. (I-18) of the dynamic 


Eq. (I-14) is a generalization of the static solution Eq. (I-12) 


of the static Eq. (I-10), incorporating time retardation. 


earch 


PUK) 


dd 
FIG. 39 Charge © moving on path r (t!). 
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From Eqs. (I-8,9,14, and 18), the integrals for b and A are 


b= | { Aut) St-te Ric) ak! at! (2-19) 


Rat) =| AG AD S(t ta R[C)AK AT (120) 


2. LIENARD-WIECHERT POTENTIALS 
The Lienard-Wiechert potentials are given in the literature 
for moving point charges of fixed magnitude. These potentials 
can be extended to the case of a fixed-position charge of time- 
varying magnitude, as shown in this appendix. 
Consider first a fixed-magnitude charge @ meving on a 
prescribed path ¥(t') with prescribed velocity (see Fig.39) 
Tit'\=c Blt’) (1-21) 
The current density at a point Y due to © is (Ref. 1, p. 465) 
F(X tVSPV sec HAYS (RS FAY] t-2D 


The vector delta function in Eq. (I-22) is defined such that 


{tex VS KS PNA = FC (1-23) 


if volume V contains point r, otherwise the integral is zero. 


Substitution of Eq. (I-22) into (I-20) gives 


(Xt) eo) | 8 5 eYS(t-teRIC)d 


(I-24) 
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¢ 
Integration of Eq. (I-24) over X; gives 


STa\- oe ( BED g(tite ReeVic) dt! (1-25) 
R (4) o| 2 8 


where 


a! eles - 
Rit = [x- eos ie 2g zi =e (1-26) 


whine 
Similarly, write for the charge density at rue) 


HRS es Cx- Fey] (1-27) 


which, when substituted into Eq. (A-19), and integrated over 
( 
Ke yields 


OK t) = e\ 1. S$ (t-tt RAVC) dt! (I-28) 


Consider now a charge in fixed position at r’, but with 

time-varying amplitude] ct’). Such a charge must, of course, 

be supplied by a current, to satisfy the continuity equation 
=-OPlot. Then the charge density at F(t!) can be written 


as 


pit) qt \$ [R- rit) | (I-29) 


Substitution of Eq. (I-29) into Eq. (I-19) gives 


Peet) =| | an) Seite Ric) § (K+ FeV a" at 


(I-30) 


The veetor potential A for q (t ) is zero, Since F(t } has zero 
velocity 

Returning to fixed-magnitude moving charges. Eqs, (1-25.28) 
are now reduced to the standard Lienard Wiechert forms 
(Ref. 1, p. 465) Write Eq. (I-25) in more convenient notation 


as 


A = el Bt) SCH N-t\ dt Gea 


where 

FEVS t+ RUtVc (1-32) 
The integral in Eq. (I-31) is evaluated using the following 
property of the delta function (Ref. 2, Pp» 465) 


{ GUNS CFG\-a)\ dX = JO. 


(I-33) 


were 


ae Fu\=a 


Hence 


= Bt d (tere ic (1-34) 
A = Ri / dt! : 


Where Eq. (1-34) must be evaluated at 


teRitVie=t 


1.€., at the retarded time 


tf ot-RtVc (I-35) 


‘See Appendix VI, 
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Afr) _ d ) ~ 1 eh ARAN (1-28) 
qe ae (Fe RAV C) = Ute <r 


The derivative ARK’ / At’ is found from 


Rit \=[ > On, - tte) ie (1-37) 


txt 


which shows RK) as an explicit function of Ke . Hence 


at ~ © axl ate 


ee we Oke 
Sa re 


om 6 
Define the unit vector AN from source charge at X to field 


point at X . ‘Then 


M\ = RN /RUEY (1-38) 
Hence 
ARID - - AU) CY (1-39) 


From Eqs. (A-21, A-36), it follows that 


LuN= Ate =e MEAN BENZ 0 (1-10) 


Hence Eq. (I-34) becomes 
pel -e Bet ) | (I-41) 
AGE ee Le 


Comparing Eqs. (I-24) and (I-28), the solution Eq. (I-41) 
to Eq. (I-24) implies that the 


Gt\- @ \ (I-42) 
ot) Ende 


Equations (I-41,42) are the standard Lienard-Wiechert potentials 


solution to Eq. (I-28) is 


which apply to a point charge of fixed magnitude © , moving at 
ty =| 
prescribed velocity w(t’) on prescribed path Y (e). 
Consider now a point charge q(t) » Fixed in position, but 


varying in magnitude with time. Integrating Eq. (I-30) first in 


Ta\ = | ged § (HVE NET reas 
BN AD 8( ‘ oud 


t 
¥ gives 


Using Eqs. (I-33, I-40), and the fact that for this case 


BUYz2TEHczo , so LUNZL 


Eq. (I-43) yields 
(f | au. (I-44) 
aad=LSe |, 


Equations (I-41, I-42, I-44) can be combined to give the 


following generalized forms of t}32 Lienard-Wiechert potentials: 
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nN f, | sinh 4 V3 (t') (I-45) 
OF RV LID 7 


gah =| =a 1-46) 
RANLEY | 


For a fixed-magnitude charge qitS= © , while for a stationary 


charge LUSN=1 » SO that Eqs. (I-45, I~46) reduce in each of 


these cases to Eqs. (I-41, I-42) and Eq. (I-44) respectively. 
Equations (I-45, I-46) were derived by Maddocks (Ref. 3). 


In MKS units, and using the notation of Ref. 2, Eqs. (41, 


42) become 
RL t\ = Aes as (1-47) 
AG) 47 { 9 
ref . 
PN enooe oe tN. 
OG, = <|s | (I-48) 
ret 
where 


ec} (I-49) 


—_ 


and ¥ is used in place of R= 


a) 
x 
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APPENDIX II 


THE FIELDS OF MOVING CHARGES 


1. EQUATIONS FOR THE FIELDS 

The fields of moving charges may be derived from the 
Lienard-Wiechert potentials Eqs. (1-41,42) and the differential 
operations of Eqs. (I-5,6) (Ref. 1, p. 466 in CGS units; Ref. 2, 
p. 345 in MKS units). The derivation below follows Ref. 1 
closely, and is included here only for the sake of completeness. 


From Eqs. (I-5,25) 


: ( 
BCH Dy =4 Ve, x Bee seer te wie /8 ]At cra 


The Vi; operator is in Xe coordinates and operates only on the 


R inside the integral. Using 
Va ea Fai\|= Veh OY KEG) cr-2) 
Eq. (II-1) becomes 


= ( 
Bus t) = e| Yi |S (t-t+ Ric VR] xBlt'\dt (II-3) 


Define 


g(R)SS(t-t eric) (R (11-4) 


- 117 - 


Using Eqs. (1-37,38) 


agce) ~ 29(8) OR a9C®) | Xi Ke! 


Dee 28|-—CU ROK aR R 
ae 
1 9(8) = >. agiey dy: = o3 & (II-5) 
Cr OXRe ‘ ro) 
Thus the UG operation is equivalent to 
n oO 
Nac= “SR (12-6) 


where mM is the unit vector along Ro.) Hence 
A t 2 
J. [SUi-te ele /R] = -W SA Ric) /R 
p A ge tr Rle) 


RC (II-7) 


( 
where C) denotes differentiation with respect to the argument 


of C) . Eq. (II-3) becomes 


Ba { Sle-te Ric) 1s (4-t+ RIG xm) at 
= = 


cR (11-8) 
Similarly, from Eqs. (1I-6,25,28,32, II-7) 


= - Oa 
eave at 


- lb = 


= ely [se-tenicr/®] at 


-~<18 = el 285 (K =teRic) St Sete Ble) A] at 


(II-3) 


A pt ._ oe (_ St) 2 S(t tric) dt’ 
C at (aa R ot 


(11-10) 


Using 
2 §(eiteRic)a~ S(t-t+elc) 
&t 


Eq. (II-10) becomes 


s = = e|-8 Stel $4 ote Ric) at! (11-11) 
Hence 
: ( 
E= e l-a (t-t+eRic)+ Sg mS (t tt 2) 


(II-12) 
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The integrals in Eqs. (II-8,12) are evaluated using the 
properties of the § ~function given by Eq. (I-33) and by Eq. 


(II-13) below (Ref. 1, p. 4). 


(hens Cra 4x 2G) as 


The first integral in Eq. (I-8) is reduced as follows: 


f= teRt'\ lc (I-14) 

dR = ( ( = ) 

ata te ge SU 11-15 
S% at = af / LU") (II-16) 


Then 


e| CBx4) § (eter At e| Bui) § (F-DAF 
R* R2L 


= 120 


Using 


CUbteRic}dt'= SC t sf IL 


the second integral in Eq. (B-8) becomes 


e | (ean) LEs8)s (f-t\af= £ £2(Fc)|, 


= £ 4, (£55 a) es 


nes 


foe 
The equation £=t means that fate S=t or that the retarded 


time th'=t- Rie . Eq. (II-8) becomes 
eri+t\2e "BAM Bx ih | (II-17) 
Bit) Rex et Gein ae 


In similaw Fashion the integrals in Eq. (II~-12) can be reduced 


by comparison with the integrals in Eq. (II-8) reduced to the 


forms appearing in Eq. (II-17). It is seen that 


oo ees a a mM B (11-18) 
E (Xt) = ea srs (#8 
R?L CL at RL ‘ 


re 
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<TR T 


Henceforth, the ee symbol witl be omitted; it being 
understood that the R.H.S. of equations for B and E are 
evaluated at retarded time < : 

To decompose the fields given by Eqs. (II-17,18) into 
"static," "velocity," and "acceleration" components, the 


differential operations in these equations must be evaluated. 


idi 2b A (RP) =~ A. aR 


as A. (-m ve 
~eER Make re 
Now 3 4 . 
QR. A NG = SH Stha,, --U (EV r1-19) 
at! — at! 2 ke) ki ~ 2 ar’ ¢ 
Using Eqs. (I-39, II-18), ry A = 
: mx (Mt) 
2 aL CACHE YoU 1 Ss (11-20) 
ost ~ CR L fe 
The derivative in Eq. (II-17) is ai 
d Ban 
at i sees ee Mat 
= & —) xi ra nS ba. »SYR-S | cr. -21) 
Equation (II-17) becomes 
ie eee ae 
Bet, 1A (Sat Saf ras] 
aes RzL cat i ten ( 


which, after some recombination, reduces to 


~ l22 - 


~ re \ 
B - @ [ aE = & ( B-\|xA (I-22) 


Equation (II-18) eae = 
A v=) 
= _—v MA t-— a. A) - 1 ate RL 

e = e| + at CL 7 
which, after considerable Cae reduces to 


= “ oo J Ad, /S- \en- 23) 
Gs e| saz aaa a CL fd (ge “Gy 


From Eq. (II-23) it is seen that 
fA = ic} 14, (& a 
MAE = e| rs Tou ae aL ) ae 
Hence 


B | =mx€\ (11-24) 
wot 


a A 
Thus, E, B, and MW form a mutually orthogonal set of vectors. 
The derivatives in Eqs. (II-22,23) are now evaluated. 


¢ x 
Denoting differentiation on t by a dot, 


A (8 = 8 _ _8 (atte) (11-25) 
dt 


RL R72 
From Eq. (I-40) - 
. A = Hy) ay (11-26) 
es —-mM-S 
Lo N28 
ae i re 


Using Eq. (I-39) 


Aen = aL tLe= RiM-8-A-B) +LCA-c 8) 


which, oy. some Beconenar we becomes 
ACRL) = _ RM -f +c p* ~CB- m (II-27) 
te 25 then i Peeuece to 
4.) = 2 BABB chS cA 
at’ RL RL RL? RL ¢rr-28) 
Similarly, Eq (I-23) becomes 


—_ p* 
ae e[ Ria SE As | - flair Ee 


(II-29) 


a = 
~ ik + ate aS Athe3)_ 2 _ bul) 
RP RL cai ce cRL 
Now - : 
— Sy SA L_ (LrdeA-8" = — 
4 
bie : ; = (II-30) 
Alsc 
m(m-S8)_ 8 A (m.s) 
CRA? CRI? cR 
= [and yA-L8- 88) 
CRL® 
= A ix((A-B x8 | a 
cRi3 


= deh - 


Substituting Eqs. (II-30,31) into Eq. (II-29), there is 


obtained (Ref. 1, p. 467) 


c - e| ar (at. 8-8") + 4. - Mx} 8x} 


(I1-32) 


From Eq. (II-2"), 
B= e| aur \_ (6 er VSx8 ta 


To convert Eqs. (II-24,32,33), which are in CGS units to MKS 


_m n(n x{CA-S ) «8})| 


(II-33) 


units, replace E by 


E Jame , e@ve/Jame, ane B oy Blemly, . 


For the sake of completeness, these equations in MKS units are 


given below as they appear in a standard text (Ref. 2, p. 356): 


B= KE (1I-34) 
Vo 
= e 7. a a, 
E-< ei = \a a2 
\ Plies Boe Na (11-35) 
ee ra [Cr- ) x 
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The following relationships hold between the symbols in Eqs. 


(II-34,35,36) taken from Ref. 2 and those used heretofore, taken 
from Ref. 1: 


Ref. 2 Ref. 1 
ie = R =RM 
i ae _ 
ale =Vic= (11-37) 
goy- H.F - R-SR= RU BeMV= RL 


2. DISCUSSION OF THE FIELDS 


The E field given by Eq. (f1-36) has three terms which 


depend differently upon the particle velocity B=V/c, For a 


Stationary charge A=S$ = O, LS \, and the field reduces to 
the ordinary static field (MKS units) 


n 


~ CN 
~ 4T7E& R* 


(II-38) 


For a charge moving at constant velocity, SH =O, and another 


term varying as \/ R? appears, which is termed the quasi-static 


or induction field. For an accelerating charge, S $0, and 


the third term in Eq. (II-32), which depends upon acceleration 


@ 
- 


B varies as ‘/® and therefore represents radiation. Similar 

remarks apply to the B fields. Equation (II-33) shows there is 
: ‘ : F 2 ; : 

an induction field varying as i/7R which depends upon velocity 


B, aua a eaasation field varying as VER which depends upon 


e 
- 


acceleration £# . 
3. RADIATION FIELDS FOR THE LOW-VELOCITY CASE 

Equations (I1I-32,33) apply to all possible velocities. The 
physics literature is replete with examples and applications to 
high velocity particles (relativistic case). However, the r. f. 
antenna engineer is concerned mostly with charges that move in 
or on conductors. While the field waveforms (also current and 
voltage waveforms) can move at essentially the speed of light c 
along antennas or transmission line conductors, the charges 
themselves move only at very low velocities compared to c. This 
is discussed in detail in Appendix IV. For such low-velocity 
particles, B=Vic is very small compared to unity, and the 


following approximations hold: 


BLL |\ 


(II-39) 


The radiation terms in Eqs. (II-32,33), denoted by subscript 


"a" for acceleration, become 


a = —£-[ tinted) (I-40) 
0. C2 R 
no art (11-41) 
(. 2 
Cc R vet 


where Q@2VU =CS = acceleration. 


For the special case where a is perpendicular to R 


("broadside" radiation), 
E - -=[4-| (II-42) 
m c* ik vot 


O.= ne | =s"| (11-43) 
(25 2 
c & tet 
In MKS units, Eqs (II-42,43) are 


C= =| 2] (11-44) 
Ge CF Ro dot 


8. => —£___[ a8 anh | (11-45) 
es Cc? ret 


which agree with Ref. 4 (pp. 376,387). 
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APPENDIX III 
WAVE PROPAGATION ON TRANSMISSION LINES 


AND LINEAR ANTENNAS 


Lt. TRAWTT TNS WAVES ON LOSSLESS TRANSMISSION LINES 

This section is intended to introduce traveling wave 
concepts. Voltage and current waves, which may be arbitrary 
functions of time, travel on the idealized, lossless, 


transmission line at a velocity given by 


ile eae mls (III-2) 
Jic 


where L is the distributed inductance in henrys per meter and 
C is the distributed capacity in farads per meter of line 
(Ref. 5, p. 12). For an air medinm, v = c = velocity of 
light = 3 x 10 m/s. The ratio of voltage to current is the 


characteristic impedance 


22) Lic ohms (III-2) 


The argument of traveling-wave functions ¢ ( Zt ) on such 
lines is of the form (tX 2/0 ). A step-function traveling- 


wave current is given by 


i@at\=1 Ul(t-2lr) (III-3) 


where U (t-a.) is the step function of unit amplitude fort7a , 
and of zero amplitude for tla. a step function at (2Z=0,t=0) 
will arrive at Z at timet=Z2lv. 

The above transmission line is a special case of the 
broader class of guided systems using the TEM mode of 
propagation (Ref. 6, p. 177). For such systems, assuming 
lossless conductors and dic’ectrics, the fields and their 
associated currents and voltages again may be arbitrary 


functions of time, and again travel undistorted at the velocity 


a\y Ge Pe (III-4) 


Jue 


2. TRAVELING WAVES WITH SINUSOIDAL TIME VARIATION ON LOSSY LINES 


For air media, v = c. 


A sinusoidally-time varying traveling-wave on a lossy line 
is of the form 
lz t\= Le ** cos(wt-82) (III-5) 
where & is the attenuation constant ( neps/m ) and ( is the 
propagation constant ( rad/m ). The propagation function & is 


C =ati A (III-6) 


and is given by 


ul 
‘= [CRrIwrY Co tyw)| . ene 


- Ju ~ 


where R,G are the series resistance and shunt conductance per 
meter, respectively. The phase velocity and characteristic 
impedance are now functions of frequency, whose equations are 
given in the literature (Ref. 5, Chap. 4). The various 
frequency components of any time waveform now travel at 
different velocities and also attenuate differently, resulting 
in distortion of any original waveform. 
3. WAVE TRAVEL OF STEP-FUNCTIONS AND PULSES ON LOSSY LINES 

Many of the characteristics of pulse waveforms traveling 
on linear antennas depend upon the propagation of such signals 
on linear transmission lines. The propagation of pulses or 
step-function wavefronts on lossy lines has been investigated 
(Ref. 5, Chap. 7). Such waveforms imply wide-band frequency 
components. There are at least two approaches to this problem. 
The first assumes that L, C, R, and G are frequency-independent. 
The second approach, which is claimed to be more realistic 
(Ref. 5, pp. 5 and 112), includes skin effects, and assumes that 
the resistance and the inductance arising from internal flux 
linkages vary as the square root of frequency, the same as the 
high-frequency asymptotes of these parameters. G is assumed to 
be ignorable, and C is assumed constant. 

For the solutions (based upon Laplace transform methods) 


and graphs of the behaviour of step-function voltages and 


currents on such lossy lines, details are available (Ref. 5, 
Chap. 7). In this appendix, some conclusions will be reviewed. 


For both of the approaches mentioned above, a time delay 
Taelay Sie (III-8) 


elapses before the voltage (or current) at distance L from an 
applied step function becomes non-zero. The velocity of 
propagation of a step~function may te said to be yas re 
only in the sense that this governs the elapsed time between 
switch-on time and the time when some effects begin to take 
place at distance L away. The calculation of these effects 
(i.e., rise time) depends upon which of the two approaches 
mentioned above is taken. 

For the variable-parameter approach, at the instanz of 
arrival of the wavefront at distance Q from the source, the 
solution predicts a smoothly changing rise in voltage and 
current, with finite rates of rise inversely proportional to 
ilies Thus the practical result of skin effect is to produce 
finite slopes of rise and fall of applied square pulses. 

4, VELOCITY OF FROPAGATION ALONG A THIN WIRE ANTENNA 

The current preduced on an antenna by a transmitter is 

a boundary value problem. Several methods have been used to 


solve this problem (Ref. 7, pp. 176-180). In this section, 


ot 


a method due to Pocklington (Ref. 8) wili be reviewed briefly, 
following Jones (Ref. 7). The details give insight into the 
phenomenon of wave propagation along a wire antenna. 
Pocklington was able to show that the current on a thin wire 
is sinusoidally distributed for sinusoidal time excitation, 
and that the velocity of the wave % the velocity of light c. 
The mathematics involved will be described mostly in words. 
The field E is written as an integral of the current, which in 
turn is used to set up an integral equation for the current. 
Assume a perfectly conducting body S (antenna less feeder aud 
receiver) in an external applied field E. @ iwt - The 
resultant total field E GNPN ae ie cengcoe Eg and an 
integral of the surface current linear density ae over . 
By using the boundary condition that the tangential component 
of E =O onS , the integral equation for a is obtained. 
Restriction to thin wires, with radius a small compared 
to antenna length and to the wavelength, reduces the integral 
equation to one dimension. The wire need not be straight, but 
sharp changes in curvature are not allowed. To find how a 
disturbance propagated along the wire, Pocklington set E,=o , 
and retained only the dominant terms in the integral equation 


for Js - To this approximation, the solution for the antenna 
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current T(L\= 214 J. at distance a) from the origin is found 


to be 
a ' yor 
Tio\= (ne 4 pele! (III-9) 
where 
raw ue (III-10) 
Using 


Ke eit pet] = co5 (wt- 2) (III-12) 


the velocity of propagation is seen to be 


wo we = oe) an (III-12) 
«K Jue 
Thus, for air, the propagation velocity along the wire is v=. 
5. SOME CHARACTERISTICS OF WAVE PROPAGATION ALONG A THIN-WIRE 
LINE 
A traveling wave of current along a thin~-wire line extend- 
ing along the Z-azis may be represented by 
Tl2t)=Tlt-2le) (III-13) 
for propagation in the positive 2 direction. The velocity of 
propagation has been taken as essentially equal to that of light 
c. The current I is related to the current density J (amps. per 
sq.m.) by 


Ti@d2| Tah 48 peiiea 
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where S is the cross-section of the wire. Determination of the 
variation of J with depth is a boundary-value problem which is 
not addressed here. The integral of J over the cross-section 
in Eq. (III-14) may be replaced by the product of a mean J taken 
as constant over some cross-section area A characteristic of the 
particular problem. Thus 

L(2,6\5 Tela, & (111-15) 
For example, for a round wire, A might be taken, as a first ap- 
proximation, as the cross-section area of the layer of thickness 
§ beneath the wire's surface where ri is a mean skin depth. Then 

TAO) = TCt-21C> crete ie) 

We now investigate the magnitude of the ratio of (POVd8 

to (WaP/dt), where is the electron charge density, and VU is 
the electron drift velocity. These quantities arise when the 
time derivative of I(Z,t) is used in calculation of radiation. 
The term (POW(dt) represents charge acceleration in the Z- 
direction and will be shown to be much greater than the other 
term (1S OHP(Ot), which represents time-rate-of-change of charge 
density. 


Using 


JFropv (ILI-17) 


there results 


— a . Sap OF ey ey (III-18) 


The continuity equation WS 2-o9lat reduces, in this one-dimen- 


sional case to 


ox .. oP (III-19) 
Oz at 

From Eq. (III-16) 
ORY - er (III~20) 
oz c 


i. : ‘ P F 
where I is the derivative of I with respect to its argument 


(t-Z/c). Likewise 


od _ Re: (TII-21) 
ot 
Combining Eqs. (III-20,21) 
-cor = ot (III-22) 
oz at 


Substituting Eqs. (III-19,22) into Eq. (III-18), 


ot ot at 


' p law /at) = £(\- = \ (III-23) 
wiapiat) wr ¢ 
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Since WlCZz| , as shown in Appendix IV, then 


plow /et) + Sf 35) (III-24) 
(arias) ~ T” 


Thus the "acceleration" texm ( PdV/dt) dominates the (woP/dt) 


term in oT fat by many orders of magnitude (about 108), That is, 


-" 9X nu 5 av (III-25) 
A oe eS at 


6. SURFACE CHARGES 

From Eq. (III-19) it is seen that there is no change of 
charge density Sowherever the current is uniform (i.e., cons- 
tant) in Z. However, if the current is not uniform along Z 
(at fixed time t), then Eq. (III-19) shows that there must be a 
change in.Pwith time. Such an antenna can be represented as a 
chain of small current elements or Hertzian dipoles each of 
length dl, with end charges which change differentially from 
element to element (Ref. 6, p. 322). These adjacent charges do 
not cancel, and the non-cancelled charges appear on the surface 
of the wire, where they produce an electric field normal to the 
wire's surface (Ref. 6, p. 322). This is illustrated in Fig.ug, 
based upon figures in Ref. 6 (p. 322). For a Hertzian dipole, 
discussed in Appendix V, the radiation fields can be obtained 
from the longitudinal (Z2-direction) acceleration of the charges 


in the wire, or equivalently, the time-derivative of the current. 
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FIG. 40 Chain of dipoles and current, charge and 
local E field distributions. 
S - Lax - 


The continuity equation automatically takes into account the net 
charges present at the dipole ends. This comment holds equally 
well for a finite-sized linear antenna with non-uniform current 
distribution, by superposition of the elementary dipoles. 

Some comments are now made about surface charges on a bent 
thin-wire line. According to Torre and Longo (Ref. 9, p. 391) 
the current and current density J must follow the bends, and 
therefore the internal E field given by J=6E must also follow 
bends. Thus electric flux follows the bends with the current 
(Ref. 9, p. 391). In order for this to be possible, additional 
surface charges must appear at the bends, as shown in Fig. 4l, 
(Ref. 9, p. 391). The net result of this, for the time-varying 
case, has not been investigated by this writer, except that it 
appears clear that the centripetal acceleration experienced by 
charges moving in the wire around bends must be due to the fields 
produced by these surface charges. These remarks apply strictly 
to a generator-driven loop, and not a loop being driven by an 
externally-applied changing magnetic flux. In the latter case 
the induced E field has nonzero curl, and forms closed loops 


without requiring any charges. 
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FIG. 4l 


Surface charges at a bend ina 
conductor carrying current. 
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APPENDIX IV 


DRIFT VELOCITY OF ELECTRONS IN CONDUCTORS 


1. FREE ELECTRON THEORY OF METALS 

According to tue free electron theory of metals, a 
conductor is a lattice of positive ions among which some 
electrons are free to move under the influence of external 
fields. In monovalent metals such as copper, siaver, and gold, 
the number of electrons averages about one per atom (Ref. 9, 
p- 144). For copper this constitutes a free electron density 
Nz iG bee va (Ref. 9, p. 147) 
2. RELAXATION TIMES 

A net (nonzero) charge density © (coulombs per mi ) 
cannot exist permanently inside a conductor (Ref. 10, p. 15, 
also Ref. 11, p. 78). In the absence of an external time- 
varying field, any pre-existing charge inside a conductor decays 
exponentially and flows to the surface. The transient decay 


equation may be derived as follows (Ref. 10, p- 15): 


J aGe (IV-1) 


Wet = of (IV-2) 


Equation (IV-1) is Ohm's law, which will be modified subse- 


quently, and Eq. (IV-2) is the continuity equation. 
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Hence 


Thus the differential equation for FP is 
OP € p=o (Iv-3) 
et «& 


The solution of Eq. (IV-3) is 
(= pldyexe.(- A=) (rven) 
at eo *\™ 8/6 


The decay time ( E/6 ) is called the relaxation time constant 
(Ref. 11, p. 78) T, , so 

Tan BLO (Iv-S) 
For copper Ty 2% lo7*> sec., which is far less than any 
characteristic times in microwaves or radio impulse systems. 
This decay process is entirely independent of an applied field, 
if present (Ref. 11, p. 78). 

There is another time constant referred to as the electron 
"nelaxation time" of a metal (Ref. 9, p. 358; Ref. 12, p. 292), 
which is the mean free time te between collisions of an average 
electron and the lattice ions, and is on the order of io ‘* secs. 
for copper (Ref. 13, p. 481). This is the time constant of the 
exponential decay of the electron drift velocity after all 


excitation is removed, as shown below. Again, this time constant 


ra ~ Jb? . 
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Thus the drift velocity Vy decays exponentially with time 
constant ts » if an applied field is removed suddenly, as 
previously stated. 


The Na which results from a sinusoidally time-varying 


EW\=E, eit i. (Ref, 13, p. 485% 


mgort Wys-e€ (rvs) 


which is obtained from Eq. (1V-12) by replacing J by \w. 
at od 


Since Wt evo w » Eq. (IV-15) becomes 


mVd_ 
te 
which is the same as Eq. (IV-8) for a constant field E. 


e \ 
Therefore, even up to frequencies on the order of \o'> Hertz, 


v.e8 (IV-16) 


Na and ) remain in phase with E , and Eqs. (1V-10,11) remain 
valid. This has been verified by experiments at K-band using 

Wl al - 

Wr Sklo tfodlsea (Ref. 13, p. 481). Thus Ohm's law JOE 


remains valid up to frequencies on the order of Lot 


Hertz, 
where © is the d.c. value. In turn, this implies that Vy 
is given correctly by Eqs. (IV-7,8). 

An alternative approach to the behaviour of conductors 


at high frequencies is through the concept of the dynamic 


conductivity (Ref. 14, p. 237). The force equation on the 
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L electron is written in the standard form 


Ave ALY -~-e Be : 
wm Sf t+bNio =-eFE (IV-17) 
at 


where 


velocity of a electron 


o <I 
a 


= a constant to account for the frictional 
resistance encountered by an electron, 
proportional to velocity. 


Ref. 14 defines 


4 
g =-2¢eV% --eSYNV; (1V-28) 


From Eqs. (IV-6,7) it is seen that 


5 = Fa-NeVy (IV~19) 
Multiplying Eq. (IV-17) by-@ , and summing over N electrons, 


wm 22 -bT= Nee (Iv-20) 
ot 


For the stationary (d.c.) case, | = constant and Eq. (IV-20) 


reduces to 


== we Ee (IV-21) 


For the stationary case, write Ohm's law as 


T=oGE (IV-22) 


fads 


ee of the conductivity. hence 


2 
6< Ne 


4 Equation (IV-20) becomes 


where 


- 4 ~~ 
bees For copper Y= 2.4 ALG : sec. % 10 \* sec. (Ref. 


be shown, 


Write Eq. (IV-24) as 


as VY PLAT ae 


Eq. (IV-22). For harmonic fields, E and J vary as 


. - LAT - 


understood that the 6© in Eq. (IV-11) is the stationary 


where 65 is defined as the stationary value (zero frequency) 


(IV-23> 


(IV-24) 


(IV-25) 


14, p. 238). 


From Eqs. (IV-11,25) it is seen that Tak, » provided it is 


° conductivity 64 . In other words, the derivation leading to 
Eq. (IV-11), (Ref. 13), assumes that G@=6@> . This assumption 


is valid for frequencies up to about 19'* Hertz, as will now 


(IV-26) 


which is recognized as the dynamic generalization of Ohm's law 


eset 


and 


Eq. (IV-26) becomes 

(jworttH T= 06 (I1V-27) 
Defining the dynamic conductivity 6 by 
J=aCGE (1V-28) 


it is seen that 


Ge (1v-29) 


which is analogous to Eq. (IV-15). For frequencies less than 
10 Me Hertz, WT4el, and OX 6 , so the dynamic Eq. (IV-26) 
reduces to the static version 
ToQE (IV-30) 
Hence Va may be calculated from 
a -=-NeVy =<QE (IV-32) 
4, ESTIMATION OF ELECTRON DRIFT VELOCITIES IN CONDUCTORS 

The drift velocity Vy will be calculated first for several 
limiting stationary (d.c.) cases, and then estimated for the 
time-varying case, 

As a first example consider a pair of parallel wires in 
air, each \ vam. in area (AWG No. 17), spaced 1 cm. center-to- 
center (Ref. 15, p. 126). Let the current be 6 amps, which is 
the maximum permissable for air-coolinrg, ‘nd the voltage be 
2,200 volts, which is almost equal to the voltage breakdown 


limit in air. The excess (surface) charpre on each wire, 


Lies + 


resulting from the capacity and impressed voltage, is 0.04 
microcoulombs per meter, which is completely insignificant 


compared to the free electron charge transported by the 


2 
> electrons Lm? » the 


25 _ 23 
number of free electrons per meter is \O- X io? =o ; 


current. For copper, with N¥ to 


Hence the conduction electron charge is 107 "x LOXKlIO> Z Le 
kilocoulombs per meter, as compared to the surface charge of 
0.04 microcoulombs per meter. The electron drift velocity 
follows from Eq. (IV-31): 


Jou —CortveatT T —- Ned 
area & 


io) tet, i Oo. gs ~4 

a Ne lo*?x1.6x to ®x soe 
It is seen that this value of Wy is negligible compared to the 
velocity of light Caz Arto F wis. 

3 a second example, consider the discussion of electric 
current and an example of the calculation of Vi by Sherwin 
(Ref. 16, p. 162-165). At any time + the free charges in all 
parts of a conducting loop have the same velocity VA, the 
Same magnitude of acceleration a, the same daldat » etc., 


that is, the charges move in unison. Then if M,= number of 
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‘ a) 


charged particies per meter whose flow is the current, and Va 
is their average velocity, then the current i = coulombs/sec 
passing a point is given by 

tomes (IV-32) 
Since WW, and © are assumed constant, then the linear charge 


density Pez gt is constant. Then 


ac ee Ne ANg = M, Ca (IV-33) 


—ae 
— 


at dt 


Q° 
d"t = me ao (IV-34) 
At? at 
Sherwin (Ref. 16) gives an example of a wire of diameter 1 mn, 
carrying current i = 100 amp., which is a large current for 
e ° 2 29 ~& 23 
this size wire. Then MM, =2WNxrarea % lO ko =lo- electrons 


per meter, and 


-~3 
a fk = So ee — T= E KO Mis 
A M, e lo*S x16 xI6? 


Sherwin (Ref. 16, p. 164-165) stat. * that 100 amp. would cause 

serious heating in this wire, - aS a practical matter 

conduction electron velocities are rarely above LO” “mis : 
Summarizing the above two examples, for currents which 


are uniformly distributed across the cross-section of a wire 
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(d.c. and low frequency for small wires), the order of 
magnitude of the electron drift velocity is not more than 
0.01 m/s. 

The high frequency or short pulse case ix now considered. 
Due to the skin effect, currents are not uniformly distributed 
across a conductor's cross-section. From Eq. (IV-7), the drift 
velocity Ny will vary with the distribution of the current 
density J across the cross-section. This also follows from 
Eq. (IV-8), since the field E varies similarly with the depth 
of penetration in the conductor. Hence calculation of Va 
requires knowledge of either J or E, per Eq. (It 31). An 
estimate of Va can be obtained by following the well-known 
assumption that a given wire at a given frequency carries its 
entire current uniformly distributed within the skin depth, 
Since the skin depth varies inversely with the square root of 
the frequency, the solution for a given exciting waveform will 
depend upon the frequency content of the waveform. Hence if 
we assume that a given total current I is confined uniformly 
with a skin depth determined by one of the higher frequency 
components of the waveform, the calculated J and Va will be 
larger than the actual average J and Va respectively. Since 


the Wa calculated in this manner is a high estimate of Vd, 
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and since the resultant Va is still many orders of magnitude 
less than c, as shown below, more accurate calculations for 
WAY » involving averaging, are not required. 

As an example consider a pulse shaped as a symmetrical 
trapezoid, with rise and fall times = t, sec. and with 
duration of the flat portion = T sec. The amplitude spectrum 
is (Ref. 17) 

\AG)|= (Tetyy|sinc CEtr) sinc $ (T+ th) | 
(IV-35) 
where sinc X= sin ( Ti% )/ (TY ). rorl =o , the shape is 
that of an isosceles triangle, with base width = 2 ty sec, 
and spectrum given by : 
( ACE)| = ty sinc (£ty) (IV-36) 
The first zero in LACE | as given by Eq. (IV-36) occurs at 
frequency £ al\{t, Hertz. For a relatively short pulse 
specified by T=O , Cy =O picoseconds, — 100,000 Mz. 
For copper, the skin depth 3 illlimeters) at frequency £ 
(Hertz) is 
SS C61 (1V-37) 
Hence for this example, SY 2]rktilotmm. Let us apply this 
to the case considered previously (Ref. 15) of a circular wire 


of diameter % 1 mm., carrying I = © amps. of current. For the 


4 
_ 
[Sal 
a) 

‘ 


short pulse, assume that this current is crowded uniformly 
into a crose~section consisting of a layer just beneath the 
wire's surface of thickness § . Without external cooling, 

this would be impossible, due to the increased heating compared 
to the already~marginal d.c. case. For the d.c. case, I is 

e iy aAé 

spread uniformly over the cross-section = W +. For the 
pulse case, I is assumed to be distributed over the layer 


cross-section = td § - It then follows that 


iP ie ee, ee | Pees 


(Na vac Rees Ric 4-§ 

where ® is the resistance per unit length. Equation (IV-38), 
as far as the resistance is concerned, is a Well-known result 
(Ref. 18, p. 176). It is based upon the assumption that the 
dynamic conductivity © = Oye » Which is valid up to f= 16 
Hertz, as discussed previously. For d = 1 mm. and $ = 2x6 “nn. , 
the ratio d/(48) is approximately = to? Using(Na)y. = 4x0 ~ 
m/S, as previously obtained, it is seen that Vg for even this 
extreme case is 0.4 m/s, which is still negligible compared to 
c= 3xXto8 m/s. 

It is concluded, therefore, that even for extremely short 
pulses of relatively large current amplitides, an upper~bound 


estimate of Na is about 1 m/s. 
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APPENDIX V 
DERIVATION OF THE RADIATION FIELDS OF AN 


ELECTRIC DIPOLE USING MOVING CHARGES 


The term "Hertzian dipole" is used to describe the follow- 
ing configuration (Ref. 19, p. 148): A charge +@ moves along 
a straight line about a stationary charge~C to form an 
electric dipole moment p (t) given by 

Pi) = @ L (v-1) 
where Q is the vector denoting the separation of the charges, 
drawn from -C to+©@ . For convenience, take Q along the 
@-axis of a spherical coordinate system (1, 8,0). The 


radiation fields of this dipole, in MKS units, derived from a 


Hertzian vector, are given by Sommerfeld as follows (Ref. 19, 


p. 150): 
Kh, — -2lNe_ P (t-rle) (v-2) 
“ance 
Ee - —a1e_ PC ric) (Vv-3) 
41 €,c7 0 


eo 
Here €) denotes differentiation with respect to time + , and 


p is the time-varying amplitude of p (t) . 
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Then, integrating over the interval ch-o2,h ) 


h 
(" _ Sae=-3=-3 | paz 
Yroz 43 : 


Multiplying by A, the desired equation is obtained, 


Y-o? 
where 
\ : Az 
Ash = rane | 2 (V-20) 


The dipole moment of the above isolated current element 


dipole is then 


ed = (EN = ( \ Tat) h et) 
Thus 
p (t- tle) = a (t- ric) h (V-22) 


and the radiation field, from Eq. (V-~2) is 


Wh, a SMO BUtevle) = AMA OE Ct-veyh 


ATiCY Ot 
a (V-23) 
For the harmonic oscillation case, it is assumed that 
TW) =F ei nw) 


~ Lod - 


From Eq. (V-23) the radiation field is 
jot -j Ke 
ee Lakes 3c (v-25) 


a 4nev 
which is a well-known result (Ref. 20, p. 498). To obtain 


this from Sommerfeld's results, note that the complex P(t) is 


= — Teh eet (V-26) 
Nee air © ° 


Since 


Toh jet | (v-27) 


———s 


Pf) = Awswt = Re | RCS 


it follows that 


AW Toh (V-28) 
°o ye 


Then the Sommerfeld resuat Eq. (V-5) reduces immediately to 
Eq. (V-25). 
Equation (V-25) will now be derived directly from che 
fields due to ~oving charges for the sinusoidally time-varying 
current element dipole of length "y described above. Consider 
a differential length az at position (@,0,2 ) along the dipole. 


The charge dq in this length is 


Aq = pfdz2 (Azarea\ (v-29) 


~ yh - 


which is taken as the term © in Eqs. (V-7,8). The accelera- 


tion of this charge is found from differentiating J] = PV ‘ 


Thus 
ea ddan tH carpe sare (V-30) 
dt! at dt! 
Then * . . 
oy XL. Laz (v-31) 
C4 = pradz — = 
PA 
Thus Eq. (V-7) becomes 
A, = ATC R - 


where R is the magnitude of the vector ® from element az 
te observation point P(TE.o). The total Hy field is found by 
integrating Eq. (V-32) over the dipole length Yi . Since LT is 


in the # direction, then for this time-harmonic case it is 


= ‘ 
T = rei a (v-33) 


and 
a aa yw (E-RIC) 
ee AM Je 2 deT - a, +4 


—-JKRR wot 
= Te "sine a ee masa 


By definition the dipole length NW is small enough so that the 


oe integral 
eee h : 
a = wee] _ S'S ane ay ed 'dz 
; oe 4c R 
oy Oo (V-35) 
g 
may be replaced by ; 
\ —JRU jut 
\ as —~ 

. &MCY (v-36) 

“ ¢@ 
Equation (V-36) is identical to Eq. (V-25). 

NY The E field of the dipole may be found from Eq. (V-36) and 


the far-field Equation (II-34), which in MKS units is 


ae E =. Mol Eo Cu x mM ) (V-37) 


%, This gives 

Ex jwoToh SINO Qe esl jut 

ATE crV (v-38) 
Equation (V-38) can also be derived from the generalized 

Lienard-Wiechert potentials (see Eqs. (1-45,46)) and the 
differential operations described in Appendix I, or from the 
integral formulations for ® and A piven by Eqs. (1-30,31). 
These derivations will not be reproduced here; for details see 


Maddocks (Ref. 3). 


APPENDIX VI 


SOME RESULTS OF DELTA FUNCTION INTEGRALS 


An integral which appears at several places in this 


report (e.g., Eqs. (102, I-33)), is of the following type: 


T= | gt) 5 (#l8)-3\ dx (VE-1) 
In Ref. 1 (p. 465) this appears as 
t= | (vI-2) 
aftlday iFoy=d 
This is not quite correct. From Ref. 1 (p. 4, example (5)), 
S§ 4D) = | (VI-3) 
\dt/d¥| | ec 20 
= the correct result is 
ng T= faa (vI-4) 
. \at/a¥ | lendea 
This can also be shown starting from the basic equation (Ref. 36, 
os p. 76) 
$ (av\= S(¥) (vI-5) 
\a\ 
F In Eq. (1-33), the variable is t! instead of x. Since df(4/)fat' 
cS 20 (see Eq. (I-40)), no error results in omitting the| |symbol. 
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Eq. (102) is of the same form as Eq. (VI-1), the vari- 


able being 2 , and g (2) =x ' bee 
£(2) = - ZCres8) /¢ 


; | d £2) |= \\-eso| 
ee Az = a 


d= -(t-rle) 


Using lz) =\ in Eq. (VI-4), Eq. (102) becomes 


y= Jo Sine Cc (VI-6) 
ATC | \- cose | 
at 


-2(rwsOV[o = - (t-te) 


—~ £42 ((-wse (VI-7) 
ee as ( » 


Since (+56) Z0 ; \ ike (05 | may be replaced by Cit ese) : 


For Of2s5h » Eqs. (VI-6,7) yield 


Hp uo ake Utt-tic) \tare (VI-8) 


40 DINO 
Niezre piuhee (VI-9) 
StS ae eG Ls) 


which is the same as Eq. (107). 


“ y6, ay = 


APPENDIX VII 
ILLUSTRATIONS OF USE OF EQUATION (74) 
FOR SOME STANDARD ANTENNAS WITH 


SINUSOIDAL TIME-VARYING EXCITATION 


l. TRAVELING-WAVE ANTENNA 
The traveling-wave current distribution on a terminated 
wire antenna of length L placed along the Z-axis from -L/2 to L/2 


is (Ref. 6, p. 342) 
TAH =T, 6 et Nien? 


Eq. (74) yields 
L(2 
jot @ 2 SIKR 


= AZ  (vrII-2) 


Hat <—ae_jwTne 
41%C 

~L/2 

The usual far-field approximation, Y272 , 


R <~ Y- 208 NEES) 


substituted into Eq. (VII-2) yields 


ul2 
2, “eRe sre ot er iM ACr asd) 
=e | dz 
4m “U2 


“WT ywt 
~ ime’ @ Sine __ sin Gn (\- 056) 


- ene \-@e 


U2 


20 
" — Tn oie yor _sime 2-2 009K (l-@56)L. 
8 ATV \-wse 
(VII-4) 
a = He \ Do [fo (VII-5) 


Eq. (VII-5) agrees with (Ref. 6, p. 342). 


2. HALF-WAVE LENGTH DIPOLE 

Consider a dipole along the Z-axis, centered at the 
origin, of total length 2h= 12. See Pip. 24, For sinu- 
soidal time-variation, the currents on the antenna are known to 


be closely approximated by (Ref. 6, p. 328) 


TV (2,t\ = Dy sin 6 (4-2) pint (270) (VII-8) 
rca, = In sin  (HtZVeNE ( 220) (vII-7) 


Then 
pu) (t= RIC) 
Ta t-R)V=Tyain ne aye! 
JER jot 
= Ty sine (nga eee 
ee, i wV 
. OT(2,t-Ric) = jeta sn RUNFAC Te (VII-8) 


ae 


Se 
a 
fr 
y 


Using Eq. VII-8, 
{RAC acy : 
Bret Ric)=jwI, sm k(hFA\é ee 


(VII-9) 


Introducing Eq. (74), and integrating, there is obtained 


es jK2@ \K2c036 
he = | | sn Wht 2) © ge fps 2\e az 
-\ 


(VII-10) 
€= —SW8__ ju Tm eee) (VII~12) 
~  &Tcr 7 
The integral in Eq. (VII-10) is (Ref. 6, p. 329) 
2 C05 ( £eas) 
Sin* 6 
Hence, using W=Wic , 
: us Jw t-te 
4, —-j tm _wolzese) oe! » (VII- 2) 


$ eur Sino 

which is a well-known result (Ref. 6, Eq. (10-59)). 

3. LINEAR ANTENNA WITH ARBITRARY CURRENT DISTRIBUTION I(Z,t) 
For sinusoidal time-varying current discributions such 


as discussed in subsectionsl and 2 above, the current can be 


represented as the product of two separate functions 


Tz, = Td et (vtt-13) 
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wee Wee . _ 1 
ty Wo 


Then 
J(wt-we) jK2oS 
ate j © 
ps a a 1,(zNe (VII-14) 


Using Eq. (74), in the far field 


: \luot-'K v\ MK Zs 
KW, 2am jee T(2)e€ Az 


4ucr antenna (VII-15) 


Eq. (VII-15) is a peneral vesult for sinusoidal time variation. 
It can also be shown that for amore general case, 
where I(Z,t) is not necessarily separable, and represents a 
general (non-sinusoidal in time) variation, that Eq. (74) gives 
the same result as obtained from the usual vector potential A 
formulation. Thus, the standard A formulation proceeds from 
(Ref. 6, p. 315) 
~ pal t 
ACE) — AL IC t-RI) av (VII-16) 
4m, R 
2 + a.) ° 
For the linear antenna \dV = JL A2@ G3. , and using kq. (VII-8) 


for the far-field, Eq. (VII-16) states that 


{ hs 
Af — A (ete b+ 2 W236) 42 as 


47 (VII-17) 


~~ 


Hence, using B=rVxE, 


ag —42_ 2 2 on [TGit-b+ See) 4,1 


(VII-19) 


The curl operation is with respect to the coordinates CY, 6, &) 
of field point P. The function in the brackets of Eq. (VII-18) 
is not a function of > » hence d1d =O . Further, since 

da= ad, se - a, sinO =, the function has no ® -component. Then 


the curl operation yields 


ak — 42! 4, [-2 d. (psu Z\--2 (me Z| 
~ 40 F 
(VII-19) 


where I has the argument arg = (2, t- tle +z! wselc) ‘ 


Retaining only the radiation term~(C\/r) , then 


t U 

dt, ——-Sie_ dz 292 t-2 + 2 ese) (vir-20) 
& 4ur or 

Using dT [dt =-(2Tlat) /ca . and YY R in the denominator, 


ditty -)\\\> rz, HeRIEVA2 ya 
47%cR at 


which is the same as Eq. (74), derived from the accelerated 


charge viewpoint. 
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